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Abstract 


We  report  for  the  first  time  thermal  non-equilibrium  (separate  electron  and  gas 
temperatures)  numerical  results  for  a  hydrazine  arcjet.  All  viscous  flow  properties  are  considered, 
assuming  laminar  axisymmetric  flow.  The  model  includes  anode  temperature  distribution,  and  the 
electrical  conductivity  is  coupled  to  the  flow  properties,  allowing  for  a  self-consistent  current 
distribution.  The  numerical  solution  algorithm  employs  the  compressible  form  of  the  PISO 
algorithm  to  solve  the  continuity  and  momentum  equations.  Run  time  is  a  few  hours  on  a  Convex 
C240  Mainframe  with  a  44  x  24  grid.  Numerical  results  are  presented  for  low  power  hydrogen 
and  hydrazine  thrusters. 

Preliminary  results  of  quadruple  electrostatic  probe  measurements  at  the  exit  plane  of  a  1 
kW  hydrazine  arcjet,  including  rie  and  Te  profiles,  are  presented.  The  quadruple  probe  model 
includes  the  effects  of  Te  and  ne  gradients  across  the  probe  volume  to  extract  Te  and  ne  radial 
profiles  from  the  asymmetric  raw  probe  data.  A  time-of-flight  electrostatic  probe  technique  for 
measuring  heavy  particle  velocities  is  described  which,  when  coupled  with  the  quadruple  probe 
data  can  yield  radial  profiles  of  ne(r),  Te(r),  Ti(r)  and  Ui(r) 

Experimental  investigations  of  the  energy  deposition  processes  in  the  nozzle  and  constrictor 
regions  of  a  1-2  kW  hydrazine  arcjet  are  being  performed.  Electron  number  density  and  electron 
temperature  measurements,  using  an  array  of  flush-mounted  Langmuir  probes,  will  be  made  in  the 
boundary  layer,  while  emission  spectroscopy  will  provide  property  measurements  in  the  core  flow. 
Current  density  measurements  will  also  be  made  with  the  pin  array,  and  will  be  compared  with 
code  results.  Application  of  a  steady-state  magnetic  field  can  be  used  to  control  arc  attachment  in  a 
predictable  manner. 
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Overview/Introduction  —  Scientific  Deficiencies  and  Barriers 

The  primary  interest  of  the  Air  Force  for  arcjet  propulsion  is  with  storable  propellants.  The 
proposed  research  emphasizes  gaseous  mixtures  of  N2  and  H2  for  the  simulation  of  N2H4  and  H2 
arcjets.  Because  the  arcjet  is  a  low  Reynolds  number  device  (Re<1000),  viscous  behavior 
dominates  energy  transfer  processes  in  the  thermal  and  viscous  nozzle  boundary  layers,  which 
couples  with  multi-dimensional  heat  conduction  in  and  from  the  massive  conducting  anode  wall. 
Therefore,  energy  transfer  must  be  described  by  coupled  Navier-Stokes  conservation  relations 
which  depend  in  a  complex  way  on  the  fluid  properties.  To  date,  no  single  computational  model 
can  be  said  to  describe  accurately  the  entire  physics  of  arcjet  flows. 

The  energy  transport  from  the  plasma  located  in  the  equilibrium  constrictor  region  of  the 
arcjet  is  a  strong  function  of  the  number  density  and  distribution  of  electrons  from  ionized  species. 
The  degree  of  non-equilibrium  can  significantly  alter  the  expected  electron  number  density.  The 
diffusion  of  species,  particularly  electrons,  is  likewise  important.  The  electric  current  distribution, 
not  known,  is  strongly  coupled  to  the  ne  and  Te  distribution  in  the  vicinity  of  the  relatively  cold 
constrictor  and  anode  walls.  Therefore  non-equilibrium  plays  a  dominant  role,  resulting  in  spatial 
variation  in  the  current  attachment  and  anode  heating.  The  proper  description  of  this  aspect  of 
arcjet  nozzle  flow  has  proven  particularly  difficult  to  achieve  in  numerical  models  of  the  arcjet. 

An  applied  magnetic  (B)  field  is  known  to  alter  electron  drift  and  hence  the  current 
distribution  and  arcjet  performance.  Because  of  these  interactions  and  because  of  the  effect  on 
swirl  in  the  constrictor,  a  numerical  model  with  applied  B  field  can  describe  the  complex  physics  of 
flow  and  energy  transfer  in  this  region  of  intense  plasma  heating.  We  expect  to  begin  research  in 
this  area  during  the  contract  year. 

All  of  the  above  factors  vary  significantly  with  power  input  (the  scaling  problem).  Note 
that  a  given  arcjet  operates  over  a  relatively  narrow  power  range,  and  DC  constricted  arcjets  as  a 
class  cannot  be  operated  efficiently  below  0.5  kW. 

This  technical  report  summarizes  the  research  carried  out  by  our  group  concerned  with 
arcjet  plasma  propulsion.  Section  I  describes  the  non-equilibrium  plasma  model  and  a  detailed 
formulation  for  kinetic  non-equilibrium.  Section  II  describes  the  experimental  research,  where 
measurements  of  electron  number  density,  electron  temperature,  Mach  number  and  gas  velocity  are 
carried  out  by  advanced  Langmuir  probes.  Section  III  outlines  our  research  with  a  modified  1-2 
kW  class  arcjet,  where  surface  Langmuir  probes  are  used  to  investigate  the  constrictor  and  nozzle 
current  attachment  regions. 
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Section  1:  Two  Temperature  Modeling 


Nomenclature 

Ci  mean  thermal  speed  (m/s) 

C  coulomb  charge 

D  effective  di^sion  coefficient  (m^/s) 
elst  elastic  energy  transfer  (W/m^) 

h  Planck's  constant  (J-s) 

hi  enthalpy  of  species  i  (J/kg) 

I  integrated  total  current  (A) 

j  current  density  (A/m^) 

kb  Boltzmann's  constant  (J/K) 

kc  Coulomb  constant  (Nm^/C^) 

mi  mass  of  species  i  (kg) 

M  Mach  number 

ni  number  density  of  ith  species  (m*^) 

p  pressure  (N/m^) 

Pq  ohmic  power  deposition  (W) 

Qij  collision  cross  section  (m2) 

Ta  anode  temperature  (K) 

Te  electron  temperature  (K) 

Tg  heavy  species  temperature  (K) 

u  axial  velocity  (m/s) 

U  heavy  species  total  particle  energy  (J) 

V  radial  velocity  (m/s) 

V  velocity  vector  (m/s) 

V  voltage  (volts) 

w  azimuthal  swirl  velocity  (m/s) 

x,r,6  cylindrical  coordinates 

Xi  mole  fraction  of  species  i 

Xmoi  mixture  parameter  in  XmolN2  +  H2 

Zi  thermodynamic  partition  functions 

ey  dissociation  energy  (J) 


*Ph.D.  Candidate,  Thomas  W.  Megli,  leads  our  group  in  developing  a  comprehensive  model  of  the  arejet 
thermophysics. 
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A,h  heavy  species  thermal  conductivity  (W/m-K) 

Ac  Spitzer  logarithm 

p  gas  viscosity  (kg/ms) 

Ei  ionization  energy  (J) 

Xa  anode  thermal  conductivity  (W/m-K) 

Xe  electron  thermal  conductivity  (W/m-K) 

V  collision  frequency  (s’*) 

p  mass  density  (kg/m^) 

a  electrical  conductivity  ^mhos/m) 

d>h  viscous  dissipation  function  (W/m^) 

Qij  collision  integral  (m^) 


1.1  Introduction 

Arcjet  thrusters  offer  a  significant  advantage  in  specific  impulse  and  therefore  cost  over 
conventional  satellite  propulsion  systems.  Lx)w  power  1-2  kW  arcjets  have  been  qualiHed  for 
north-south  station  keeping,  while  higher  power  20-30  kW  designs  are  being  considered  for  orbit 
transfer  applications.* 

Most  arcjet  design  strategies  to  date  have  been  empirical.  The  need  for  a  more  complete 
description  of  the  complex  physical  processes  in  an  arcjet  has  resulted  in  a  variety  of  numerical 
models  ranging  from  simplified  one-dimensional  design  tools^  to  more  complex  two-dimensional 
descriptions  of  the  gas  dynamics  and  energy  transfer  processes.3-9  Comprehensive  arcjet  models 
are  required  to  (1)  interpret  experimental  results,  (2)  understand  the  physical  processes  in  regions 
of  the  thruster  where  diagnostic  techniques  are  challenging,  (3)  improve  arcjet  performance  and  (4) 
ultimately  scale  the  devices  to  higher  power  levels. 

The  fluid  dynamics  and  energy  transfer  within  the  arcjet  thruster  are  nonlinear  and  strongly 
coupled.  The  energy  conversion  mechanism  for  an  arcjet  is  an  electrical  discharge,  or  arc, 
produced  by  a  voltage  difference  between  the  thruster  nozzle,  which  is  the  anode,  and  a  conical 
electrode  on  the  upstream  side  of  the  constrictor,  which  is  the  cathode.  The  arc  current 
distribution  depends  on  several  factors  in  addition  to  the  geometry.  The  distribution  is  coupled  to 
the  propellant  mass  flow  rate,  composition,  thermal  properties,  electrical  conductivity,  and  gas- 
dynamic  properties. 

Both  chemical  and  thermal  processes  in  the  arcjet  are  described  as  non-equilibrium.  Arc  current 
is  converted  to  electron  thermal  energy  through  ohmic  dissipation.  The  electrons  transfer  thermal 
energy  to  heavy  species  in  the  arc  plasma  through  collisions.  This  energy  is  then  converted  to 
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thrust  as  the  fluid  accelerates  through  the  nozzle.  In  regions  of  low  pressure  and/or  elevated 
ohmic  heating,  collisional  coupling  between  electrons  and  heavy  species  may  not  establish  equal 
gas  and  electron  temperatures.  Additionally,  flow  velocities  are  large  so  that  fluid  residence  times 
in  the  nozzle  are  of  the  order  of  1  ps;  much  of  the  energy  invested  in  dissociation  and  ionization  is 
frozen.  Pressures  drop  from  roughly  1  atm  at  the  constrictor  to  millitorr  conditions  at  the  exit; 
radial  diffusion  of  electrons  from  the  arc  core  is  important  in  determining  the  arc  structure. 

The  thrust  and  specific  impulse  produced  by  an  arejet  are  determined  by  several  factors  including 
the  power  transferred  to  the  propellant,  the  extent  of  both  kinetic  and  chemical  non-equilibrium, 
the  propellant  gas  mixture,  and  the  arejet  nozzle  geometry.  Other  factors  which  contribute  to  the 
overall  performance  include  heat  loss  to  the  anode  and  the  voltage  sheaths  at  the  electrodes.  A 
comprehensive  model  of  an  arejet  thruster  is  therefore  highly  complex,  including  plasmadynamic, 
fluid  dynamic,  and  surface  and  volumetric  heat  transfer  phenomena. 

This  paper  describes  the  formulation  of  a  self-consistent  arejet  model.  In  section  1 .2  the  kinetic 
non-equilibrium  model  is  described,  and  in  section  1.3  the  numerical  solution  method  is  presented. 
Numerical  results  are  presented  in  section  1.4.  Non-LTE  solutions  are  presented  for  both 
hydrogen  and  hydrazine  cases. 

The  complete  arejet  model  will  ultimately  include  detailed  descriptions  of  kinetic  non-equilibrium, 
finite  rate  chemistry,  radiation  transfer,  and  species  diffusion.  The  model  discussed  here  focuses 
on  the  kinetic  non-equilibrium  aspects  of  arejet  thrusters,  and  will  provide  a  baseline  comparison 
for  future  research  including  both  kinetic  and  chemical  non-equjlibrium. 

1.2  Model 

An  axisymmetric,  steady  state,  laminar,  two  temperature  kinetic  non-equilibrium  model  is 
formulated  for  a  direct  current  arejet  with  flow  swirl,  a  variable  nozzle  geometry,  and  a  variable 
mixture  of  nitrogen  and  hydrogen.  The  model  predicts  thrust,  specific  impulse,  and  internal  fields 
for  pressure  p,  density  p,  and  both  electron  Te  and  heavy  species  Tg  temperatures.  Additionally, 
the  model  also  predicts  the  cuirent  density  distribution  j,  voltage  potential  V,  and  velocity 
components  u,  v,  and  w. 

The  details  of  the  model  are  presented  in  six  sections  below  as:  (1)  the  fluid  dynamic  equations, 
(2)  the  energy  equations,  (3)  the  electromagnetic  equations,  (4)  the  species  number  density 
computation,  (5)  transport  coefficients,  and  (6)  boundary  conditions.  Since  an  azimuthal  velocity 
component  w,  or  flow  swirl,  is  included  in  the  model,  an  azimuthal  momentum  equation  must  be 
included  and  a  term  containing  w  is  included  in  the  radial  momentum  equation.  The  arejet 
geometry  is  described  as  (x,r,0)  where  x  is  the  axial  coordinate,  r  is  the  radial  coordinate,  and  0  is 
the  azimuthal  coordinate. 
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The  viscous  fluid  dynamic  equauuns  are  summarized  below: 


Axial  momentum: 


a^puu)^j^a(rpuv)  dj  + 

dx  r  Sr  3  dx\  dx|  r  9i^  3r|  3  dx\  dr)  3r  dx  r  dr\^dx|  dx 


Radial  momentum: 


9(puv)  I  a(rpvv)  _  ^  2v?E 

Bx  f  ar  3x1  3x1  3r  drl  3r  /  3^  Br }  3r  Bx  3r  ar 

3  r2  3r  dn  ox)  ar  r 


Azimuthal  momentum  (Bp/BO  s  0): 


Axisymmetric  continuity: 


|^pu)  +  7|(rpv)  =  0 


The  assumption  of  kinetic  non-equilibrium  gives  separate  energy  equations  for  the  electrons  and 
heavy  species.  The  electron  energy  equation  is 


r  •  V  •  • 

V  •  ^kb(neV  -  |^)Te  =  V  (XeVTe)  +  iiieV  (heDeVne)  -  PeV  v  +— -  -  elst  -  radiation  (5) 
V2  c  /  G 


where  'elst'  represents  energy  lost  through  elastic  collisions  with  the  heavy  species,  and 
'radiation'  denotes  optically  thin  radiation  loss  due  to  continuum  bremsstrahlung  radiation.^  ^  The 
elastic  energy  transfer  source  term  is  calculated  as  the  sum  of  electron-atom  and  electron  ion 
contributions  as 


elst  =  3kb(Te  -  Tg)nemc  I 


heavies 


where  the  average  collision  frequencies  between  electrons  and  heavy  species  are  calculated  using 
the  mean  electron  thermal  speed  and  collision  cross  s^tions  Qes  as^^ 


(7) 


_  /SkbTe  „  ^ 

.  ^sQes 
7Cme 


Ves  =  '/ 


where  the  electron  ion  cross  sections  are 


^  _  Tce^k^lnAc 
vei  - 

2(kbTc)2 


(8) 


Electron-H2  and  electron-N2  collision  frequencies  are  multiplied  by  energy  loss  factors  to  account 
for  inelastic  rotational  and  vibrational  losses.  •  7 
The  heavy  species  energy  equation  is 

V  (nhvU)  =  V  (XhVTg)  +  2,  mjV  (hiDiVni)  -  phV  v  +  elst  +  Oh  (9) 

i  ^  e 


where  the  particle  energy  U  includes  translational,  rotational,  vibrational,  and  bulk  kinetic  energy 
contributions  as  well  as  ground  state  ionization  and  dissociation  potentials. 

The  anode  temperature  distribution  is  determined  from  the  energy  equation  for  heat  conduction  in 
a  solid.  Neglecting  ohmic  dissipation  in  the  anode  gives 

v(XaVTa)  =  0  (10) 


where  X-a  =  MTa)^^. 


3.  Species  Equations 

The  species  number  densities  are  computed  iissuming  a  tv/o-temperature  chemical  equilibrium. 
With  this  formulation  the  species  number  densities  are  given  as  a  function  of  the  pressure  p,  the 
heavy  species  temperature  Tg,  and  the  electron  temperature  Tg.  A  plasma  composition  of  N2,  H2, 
N,  H,  N+,  H+,  and  electrons  is  assumed.  Dissociation  and  ionization  equations  are  constructed 
from  a  generalized  law  of  mass  action  as  given  by  entropy  maximization.  This  yields  the 
following  equations  for  the  nitrogen  and  hydrogen  species: 


'^H2 


/7imHkbTg\3/2 

f  Z^,ex  1 

A  kbTg 

'  h2  ) 

\ZH2,exZH2,rotZH2,vib  / 

(11a) 
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*^N2  \  ^bTg  ^  |j2  /  \ZM2,exZN2  ,rotZN2,vib 


(11b) 


neOH-K 

"H 


Dn 


(11c) 


(lid) 


In  the  above  expressions,  the  excitation  partition  functions  Zex  are  functions  of  Tg,  while  Zrot 
and  Zvib  are  functions  of  Tg.  Details  for  the  calculation  of  the  various  partition  functions  can  be 
found  in  References  11  and  15-17. 

The  plasma  equation  of  state,  quasi-neutrality,  and  the  definition  of  the  mixture  parameter  Xmoi 
provide  the  additional  relations  required  to  solve  for  the  seven  unknown  species.  Note  that  the 
initial  mixture  composition  is  written  as  XmoiN2  +  H2  so  that  pure  hydrogen  (xmoi  =0),  simulated 
hydrazine  (Xmoi  -  1/2)  and  simulated  ammonia  (Xnioi  =  1/3)  can  be  easily  investigated: 


2nN2  +  n^  ■^nN^. 

2nH2  +  uh  +nH+ 


(12) 


ric  =  nH+  +  nN+ 


(13) 


p  =  nekbTe  +  (nN2  +  nH2  +  )l^Tg  (14) 

Equations  (1 1)  -  (14)  are  iteratively  solved  using  the  Newton-Raphson  method.^*  Results  are 
indicated  in  Figures  1-1  a  through  1-lc,  where  the  thermal  equilibrium  mole  fractions  for  hydrogen 
(xmoi  =  0),  hydrazine  (Xmoi  =  0.5),  and  nitrogen  (Xmol  -->  plasmas  are  shown.  Hydrazine 
results  at  p  =  1  atm  indicate  nearly  complete  dissociation  of  H2  at  gas  temperatures  above  Tg  = 
5000  K,  followed  by  complete  N2  dissociation  above  Tg  =  8000  K.  The  onset  of  ionization 
begins  at  approximately  Tg  =  10000  K.  Complete  ionization  is  noted  for  temperatures  in  excess  of 
Tg  =20000  K.  These  calculations  agree  well  with  those  of  Capecchi  and  d’Agostino.l^ 
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Figure  1-la:  Hydrogen  mole  fractions  at  p  =  1  atm. 
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Figure  l-lb:  Hydrazine  mole  fractions  at  p  =*  1  atm 


p=  1.0  atm 


4.  Electromagnetic  Equations 

The  calculation  of  the  electrical  current  distribution  is  critical  to  understanding  arc  physics  and 
poses  a  challenging  problem  in  arcjet  analysis.  The  current  distribution  is  largely  determined  by 
the  flow-field  electrical  conductivity,  which  is  a  strong  function  of  the  electron  number  density  and 
temperature  distribution.  Local  thermodynamic  equilibrium  (LTE)  models  must  incorporate 
minimum  electrical  conductivity  "floors"  in  the  cool  boundary  layer  region  to  avoid  excessive 
voltage  drops  near  the  anode.^-®  The  artificial  restriction  that  Te  =  Tg  results  in  an  underprediction 
in  the  ionization  (and  therefore  electron  number  density  and  o)  in  the  boundary  layer  regions  and  a 
subsequent  overprediction  in  power  deposition.  Hence,  the  conductivity  floor  is  employed  to 
provide  agreement  between  the  experimentally-observed  and  predicted  power  deposition. 

Diffusion  of  electrons  from  the  arc  core  towards  the  anode  wall  also  elevates  the  electron  number 
density  and  plays  an  important  role  in  determining  the  boundary  layer  electrical  conductivity.  This 
has  been  demonstrated  by  Butler,  Kull,  and  King,  who  show  that  the  incorporation  of  diffusion 
provides  a  reasonable  current  distribution  without  Imposing  artificial  restrictions  on  a  or  the 
current  attachment  region.^  However,  the  two-temperature  solutions  of  Miller  and  Martinez- 
Sanchez^  demonstrate  that  significant  kinetic  non-equilibrium,  with  elevated  Te,  exists  outside  of 
the  arc  core  in  the  constrictor  region.  These  results  indicate  that  both  diffusion  and  increased 
ionization  due  to  elevated  Te  are  important  factors  which  determine  the  boundary  layer  e  jctrical 
conductivity. 

The  correct  approach  to  the  arc  attachment  problem  is  to  model  the  anode  as  an  equipotential 
surface  and  allow  the  current  distribution  to  be  a  model  output  which  is  independent  of  artificial 
restrictions.  We  are  adopting  a  self-consistent  approach  for  solving  the  current  distribution  and 
electric  field  simultaneously,  for  which  the  current  distribution  is  properly  coupled  to  the  flow- 
field. 

Assuming  induced  magnetic  fields  are  negligible,  and  neglecting  electron  pressure  gradient 
driving  terms,  ^  2  Ohm's  law  and  current  conservation  are  used  to  solve  for  the  current  density  j  and 
voltage  V: 


II 
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(15) 
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5.  Transport  Coefficients 

Transport  coefficients  for  the  present  research  are  approximated  using  mean  free  path  mixture 
rules  taken  from  Mitchner  and  ICruger.^^  The  viscosity  is  computed  as 


N 

4=1 

i=I 


N 

I 

J=1 


-4i 


njMij 


where  Mjj  is  defined  as 


-^22) 


and  the  pure  species  viscosity  is 


Pi  =  l/2mjniCili 


(17) 


(18) 


(19) 


where  cj  is  the  mean  thermal  speed,  1}  is  the  mean  free  path,  and  flijj  is  the  collision  integral  for 
interaction  between  species  i  and  j. 

The  translational  contribution  to  the  heavy  species  thermal  conductivity  is  calculated  using  an 
equation  similar  to  (17)  with  the  pure  species  conductivity  given  as: 


Xi  = 


(20) 


The  election  thermal  conductivity  is 


X^  =  .  2A. 

|l  + 

\  Ve^,/ 


(21) 


where  the  electron-ion  (e-i)  and  electron-heavy  (e-h)  effective  collision  frequencies  are  summed 
over  the  ionic  and  heavy  species. 

The  electrical  conductivity  is  calculated  from 
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me  I  niCeilc-i 

i  ^  e 


The  effective  diffusion  coefficients  Dj  are  calculated  from 


where  the  Dij  are  the  binary  mixture  diffusion  coefficients  given  as^l 


/2kbT(mi+  mj)U/2 


IbnQ 


Jcmimj 


The  electron  and  ion  diffusion  coefficients  are  modified  to  account  for  the  ambipolar  electric  field. 


Assuming 


"H+  “N+ 


,  the  ambipolar  diffusion  coeffiecients  for  electrons,  H+,  and  N+ 


***  +  ^H+nH+  +  HN+nN+ 

_  PH+(Mcne  l^N-fPN-^)  +  ltH-»(Pene  •  DN^^nN^^) 
M«ne  +  liH+nH+  +  liN+nN+ 

JJ  _  DN->(Metle  -*•  llH-^nH^)  +  HN-»(PePe  •  DH-tJH-t-) 
Me^e  +  llH+nH+  +  llN+nN+ 


where  the  mobilities  p  are  calculated  from  pi  s 


W 


The  collision  integrals  and  cross-sections  required  for  the  above  relations  are  obtained  from 
References  22-27. 

Equilibrium  transport  coefficients  for  hydrazine  are  shown  in  Figures  1-2  through  1-5.  The 
heavy  species  viscosity,  indicated  in  Fig.  1-2,  varies  widely  due  to  the  changing  gas  composition, 
with  maxima  at  Tg  =  BOOO,  10000,  and  12000  K  for  pressures  of  p  =  0.01,  0.1,  and  1  atm.  The 
maxima  increase  with  pressure,  indicating  peak  values  of  roughly  p  =  1.7  x  10*^,  2.0  x  10^,  and 
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2.4  X  10'^  N-s/m^  respectively.  Qualitatively  similar  results  for  the  translational  thermal 
conductivity  are  shown  in  Fig.  1-3. 

The  electron  thermal  conductivity,  shown  in  Fig.  1-4,  essentially  vanishes  for  temperatures  below 
6000  K.  For  T  >  6000  K  Xe  increases  monotanically.  Electrical  conductivities  are  shown  in  Fig. 
1-5.  For  5000  K  <  T  <  9000  K,  O  varies  inversely  with  pressure.  Above  9000  K,  o  increases 
with  pressure  and  temperature,  indicating  values  at  T  =  20000  K  of  O  =  10000,  12500,  and  15500 
mhos/m  for  pressures  of  p  =  0.01,  0.1,  and  1  atm. 


Figure  1-2:  Viscosity  versus  temperature  for  hydrazine. 


Figure  1-3:  Translational  thermal  conductivity. 
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Figure  1-4:  Electron  thermal  conductivity. 


Figure  1-5:  Electrical  conductivity. 
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6.  Boundary  Conditions 

The  physical  boundary  conditions  for  the  continuity,  momentum,  energy,  and  electromagnetic 
equations  are  summarized  in  Table  I  and  Fig.  1-6.  At  the  inflow  boundary  lA,  the  flow  is  axial 
and  the  stagnation  temperature  Tt  and  stagnation  pressure  pt  are  specified.  The  stagnation  pressure 
is  adjusted  to  obtain  the  desired  mass  flow  rate.  The  static  pressure  p  is  extrapolated  from  the 
interior  of  the  domain.  The  u-velocities  and  the  inlet  static  temperatures  are  then  computed  from  the 
isentropic  stagnation  property  relations.  The  effects  of  flow  swirl  are  incorporated  by  specifying 
the  azimuthal  inlet  w-velocities  as  a  fraction  of  the  u-velocities.  At  the  solid  cathode  1C  and  anode 
ID  surfaces,  no-slip  conditions  are  employed.  Zero  radial  gradients  are  specified  on  the  center-line 
IE. 

At  the  outflow  boundary  IB,  since  the  boundary  layer  flow  is  subsonic,  a  low  static  pressure  is 
specified.  This  is  incorrect  for  supersonic  flow.  The  effect  of  this  condition  on  the  solution  will 
be  systematically  investigated;  the  correct  exit  pressure  distribution  would  ultimately  require  the 
extension  of  the  computational  domain  to  include  the  plume.  We  are  also  considering  a  mixed 
boundary  condition  where  the  pressure  is  extrapolated  for  locally  supersonic  flow,  and  is  fixed  for 
locally  subsonic  flow. 

The  outflow  u-velocities  are  updated  to  satisfy  continuity  using  velocity  corrections  as  given  by 
discretized  forms  of  the  momentum  and  continuity  equations  for  the  exit  plane.  This  is  consistent 
with  the  numerical  solution  method  discussed  in  section  1.3.  Note  that  a  zero  electron  temperature 
gradient  is  employed  at  the  walls  1C  and  ID.  At  low  temperatures,  the  electron  thermal 
conductivity  approaches  zero  so  that  this  boundary  condition  will  have  little  effect  on  the  results.^* 

The  boundary  conditions  for  the  voltage  are  shown  in  Table  I.  The  anode  is  modeled  as  an 
equipotential  surface,  while  the  cathode  voltage  is  either  equipotential  or  has  a  prescribed  axial 
current  distribution  at  the  tip.  Note  that  in  this  formulation  the  sheath  voltage  drops  Vs  have  been 
uncoupled  from  the  bulk  plasma.  The  axial  gradient  of  the  voltage  is  set  to  zero  at  the  inflow  and 
outflow  boundaries  so  that  current  is  forced  to  attach  within  the  arcjet  nozzle.  Recent  radial  electric 
field  and  a  measurements  along  the  exit  plane,  showing  the  existence  of  roughly  1  A/cm^  current 
density,  suggest  that  this  boundary  condition  may  need  refinement.29 
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Table  I:  Boundary  conditions. 


Figure  1-6:  Schematic  of  fluid  and  anode  domain  boundaries. 


The  boundary  conditions  for  the  anode  energy  equation  are  also  shown  in  Table  I.  On  the  outer 
surfaces  2B  and  2D,  the  local  heat  flux  is  given  by  radiation  to  the  surroundings  at  a  specified 
background  temperature  Tbg  «  Ta.  Along  the  anode  surface  2C,  the  net  heat  flux  is  calculated  as 
the  sum  of  the  heat  flux  due  to  conduction  from  the  flow-field  and  the  sheath  losses  given  by  the 
product  of  the  local  current  density  and  the  estimated  sheath  voltage  drop  Vs.  Note  that  the 
radiation  loss  from  the  plasma  to  the  anode  is  neglected  in  this  formulation.  Surface  2A  is  either  a 
fixed  temperature  or  zero  heat  flux  boundary. 
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1.3  Numerical  Solution  Method 


The  solution  method  is  structured  around  the  global  continuity  and  momentum  equations.  These 
are  solved  using  a  compressible  form  of  the  pressure-based  PISO  algorithm,^^  where  the  density 
variations  are  implicitly  included  in  the  pressure-correction  procedure  as  described  by  Rhie.^  ^ 
Staggered  velocity  and  scalar  grids  yield  a  well  connected  pressure  field.  Upstream  weighted 
densities  are  employed  to  insure  numerical  stability.^^  The  equations  are  discretized  and  power 
law  differencing  of  convective  and  diffusive  portions  of  the  linearized  coefficients  preserves 
numerical  stability  in  regions  of  the  flow  where  convection  dominates.  The  solution  process  is 
iterative  and  values  for  the  field  variables  from  the  latest  iteration  are  used  to  linearize  the  source 
terms.  Under-relaxation  is  employed  so  that  updated  solutions  are  fractionally  added  to  the 
solution  from  the  previous  iteration.  The  governing  equations  are  transformed  into  natural 
coordinates  and  solved  on  a  uniform  computational  mesh. 

Fig.  1-7  illustrates  the  solution  algorithm.  The  energy  equations  are  solved  to  update  gas, 
electron,  and  anode  temperature  fields.  The  species  number  densities,  transport  coefficients,  and 
source  terms  are  then  updated  using  the  latest  available  temperatures  and  pressures.  The  continuity 
and  momentum  equations  are  solved  to  update  the  pressure  and  velocity  fields.  Finally,  the  voltage 
and  current  density  distributions  are  updated.  Iterations  are  continued  until  the  discretized  equation 
residuals  are  converged. 
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Figure  1-7:  Numerical  solution  algorithm. 


1.4  Model  Results 

Initial  results  for  the  arcjet  model  were  summarized  by  Megli,  et  al.^^  An  LTE  solution  was 
given  for  a  low-power  hydrogen  arcjet  flow.  The  vanishing  near-anode  electrical  conductivity  a 
inherent  in  the  assumption  of  thermal  and  chemical  equilibrium  required  the  use  of  an  analytical 
electrical  conductivity  to  avoid  unreasonably  high  ohmic  power  depositions.  Both  the  anode  heat 
transfer  and  the  current  distribution  were  decoupled  from  the  flow. 
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Recent  modeling  progress  has  yielded  thermal  non-equilibrium  solutions  including  coupled 
electrical  conductivity  and  anode  heat  transfer.  Te  >  Tg  predicts  a  finite  near-anode  electrical 
conductivity.  This  is  clearly  illustrated  in  Fig.  1-8,  where  the  electrical  conductivity  for  a  hydrogen 
plasma  at  p  =  1  atm  is  shown  for  various  combinations  of  Tg  and  Te-  With  Tc  =  Tg,  a  is 
essentially  zero  for  Tg  <  4000  K.  As  the  ratio  of  gas  to  electron  temperature  is  increased,  the 
increased  ionization  elevates  a  from  the  value  predicted  assuming  LTE  (Te/Tg  =  1.0).  For 
example,  for  Tg  =  6000  K,  a  increases  from  roughly  a  =  50  mhos/m  to  a  =  5000  mhos/m  as  Te  is 
increased  to  12000  K  (Te/Tg  =  2.0). 


Hydrogen  Plasma  at  P  =  1  atm. 


Figure  1-8;  The  effect  of  kinetic  non-equilibrium  on  the  electrical  conductivity. 

The  numerical  parameters  and  results  for  three  non-LTE  solutions  are  summarized  in  Table  II. 
For  these  cases,  flow  swirl,  species  diffusion  contributions  to  the  heat  flux  vector,  and  viscous 
dissipation  Oh  are  neglected.  The  exit  pressure  is  constant  at  p  =  25  Pa.  The  temperature  of  the 
upstream  anode  surface  2 A  is  fixed  at  T  =  1000  K,  and  the  cathode  temperature  increases  linearly 
from  T  =  10(X)  K  at  the  inlet  to  T  =  3000  K  at  the  tip.  Also,  a  constant  anode  fall  of  Vg  =  7.0  volts 
is  assumed,  and  a  uniform  axial  current  is  specified  on  the  cathode  tip.  The  emissivity  of  the  anode 
is  e  =  0.32,  and  the  stagnation  temperature  of  the  inflow  is  Tt  =  1(X)0  K. 


case 

pjt  at  inlet  (Pa) 

rh  (mg/s) 

Pn(W) 

1 

Hi 

240000 

15 

17.1 

111 

912 

2 

Hi 

240000 

20 

12.6 

1044 

1011 

3 

N2H4 

475000 

15 

48.5 

949 

613 

Table  H:  Parameters  for  non-LTE  simulations. 
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For  each  case,  a  grid  with  44  uniformly  spaced  axial  nodes  is  used.  The  fluid  and  anode 
domains  contain  11  and  12  radial  nodes  as  indicated  in  Fig.  1-9.  These  grids  clearly  warrant 
refinement  in  the  constrictor  region;  however,  the  coarser  grids  minimize  the  computational  effort 
required,  and  the  solutions  reveal  interesting  non-equilibrium  physics  present  in  arcjet  flows. 
Additional  cases  on  refined  grids  are  in  progress.  Computations  are  performed  on  the  UIUC 
Convex  C-240  mainframe  computer. 


Gas  temperature  contours  are  shown  in  Figures  1-lOa  through  1-lOc.  For  the  low  power 
hydrogen  case  1  =  45  MJ/kg),  a  peak  temperature  of  roughly  Tg  =  7500  K  is  noted  near  the 

cathode  tip.  The  gas  temperature  then  decreases  to  Tg  =  1500  K  as  the  flow  accelerates  along  the 

centerline  to  the  nozzle  exit.  The  higher  =  83  MJ/kg  indicated  for  case  2  gives  a  higher 

th 

temperature  of  Tg  =  9000  K  in  the  constrictor  and  a  higher  centerline  exit  temperature  of  Tg  =  2800 

K.  For  the  hydrazine  case  3  =  20  MJ/kg),  the  gas  temperature  in  the  constrictor  is  roughly 

m 

Tg  =  12000  K  and  decreases  to  roughly  2900  K  at  the  exit. 


Gas  temperature  contours  [K] 


^  X  [mm] 


Figure  1-lOa:  Gas  temperature  contours  for  case  1  (H2  at  =  45  MJ/kg). 


Gas  temperature  contours  [K] 
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Figure  1-lOb:  Gas  temperature  contours  for  case  2  (H2  at  ^  =  83  MJ/kg). 


Gas  temperature  contours  [K] 


X  [mm] 


Figure  1-lOc:  Gas  temperature  contours  for  case  3  (N2H4  at  =  20  MJ/kg). 
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Thermal  non-equilibrium  is  indicated  in  Figures  1-1  la  through  1-1  Ic,  where  contours  of  Te/Tg 

are  shown.  A  comparison  of  Figures  1-1  la  and  1-1  lb  indicates  that  increasing  the  power 

deposition  from  --p-  =  45  MJ/kg  to  =  83  MJ/kg  drives  the  flow  towards  thermal  equilibrium. 
'  m  m 

Increased  ionization  improves  the  collisional  coupling  between  the  electrons  and  the  gas.  The  flow 
near  the  anode  is  highly  non-equilibrium  with  maxima  just  downstream  of  the  constrictor  of 
roughly  Tg/Tg  =  18,  9,  and  8  indicated  for  cases  1-3  respectively.  The  elevated  Tg  in  this  region 
produces  the  ionization  necessary  for  a  reasonable  near-anode  a  prediction. 

Axial  velocity  contours  are  shown  in  Figures  1-1 2a  through  l-12c.  For  the  hydrogen  flows, 
peak  centerline  velocities  of  roughly  u  =  16000  m/s  and  u  =  18000  m/s  are  indicated  at  an  axial 
location  of  x  =  3  -  4  mm  downstream  of  the  cathode  tip.  Interestingly,  the  centerline  velocity 
increases  monotonically  toward  the  nozzle  exit  for  the  hydrazine  case. 

Anode  Temperature  Distributions 

Figures  1-1 3a  through  1-1 3c  illustrate  the  anode  temperature  distributions.  Increasing  the  power 

deposition  and  current  increases  the  thermal  loading,  as  peak  temperatures  near  the  constrictor  of 

roughly  Ta  =  1 190  and  Ta  =  1460  K  are  noted  for  cases  1  and  2.  This  is  clearly  illustrated  in  Fig 

1-14,  where  the  conductive  and  sheath  loss  contributions  to  the  anode  heat  flux  are  shown. 

Maxima  in  conduction  losses  of  q  =  2500  and  q  =  8500  W/cm^  occur  at  the  constrictor  entrance  for 

^  =  45  IvU/kg  and  =  83  MJ/kg  respectively.  Peak  sheath  losses  of  q  =  3000  and  q  =  6400 
m  m 

W/cm2  are  noted.  Since  a  constant  anode  fall  of  Vg  =  7.0  Volts  is  assumed,  the  sheath  losses  are 
simply  proportional  to  the  local  current  density,  and  the  maxima  occur  just  downstream  of  the 
constrictor  where  the  majority  of  the  current  attaches. 
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Figure  1-1  la:  Te/Tg  contours  for  case  1  (H2  at  =  45  MJ/kg). 
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Rgure  l-12a: 


Axial  velocity  u  contours  for  case  1  (H2  at 


^  =  45MJ/kg). 
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Figure  l-12b:  Axial  velocity  u  contours  for  case  2  (H2  at 


^=83MJ/kg). 
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Figure  l-12c:  Axial  velocity  u  contours  for  case  3  (N2H4  at 
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Figure  I- 13a:  Anode  temperature  distribution  for  case  1  (H2  at 


^  =  45MJ/kg). 
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Figure  1-14:  Heat  losses  to  the  anode  for  cases  1  and  2. 


1.5  Conclusions 

A  two-temperature  chemical  equilibrium  model  has  been  formulated  for  N2/H2  arcjets.  A 
modified  set  of  Navier-Stokes  equations  employs  separate  electron  and  heavy  species  energy 
equations  to  account  for  separate  heavy  and  electron  temperatures.  The  assumption  of  chemical 
equilibrium  allows  for  a  tractable  solution  for  the  species  number  densities  using  the  laws  of  mass 
action,  the  plasma  equation  of  state,  quasi-neutrality,  and  species  conservation.  In  this 
approximation,  finite-rate  chemistry  and  diffusion  are  neglected  and  the  species  number  densities 
are  determined  from  the  pressure,  and  electron  and  heavy  temperatures. 

Initial  non-LTE  results  indicate  that  thermal  non-equilibrium  throughout  the  arcjet  boundary  layer 
significantly  reduces  the  energy  transfer  from  the  electrons  to  the  heavy  species,  and  has  important 
implications  for  the  prediction  of  electrical  conductivity.  Increasing  the  specific  power  deposition 
reduces  the  non-equilibrium  and  improves  Isp  but  increases  the  anode  thermal  loading.  The  model 
is  capable  of  producing  results  for  both  H2  and  N2H4  arcjets.  This  work  represents  progress 
toward  a  systematic  evaluation  of  the  effects  of  kinetic  non-equilibrium  on  arcjet  performance  with 
the  goal  of  providing  a  useful  numerical  tool  for  interpreting  experimental  results. 
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Section  2:  Exit  Plane  Electrostatic  Probe  Studies* 

2.1  Quadruple  Probe  Experiments 
Introduction 

The  use  of  electrostatic  probes  has  long  been  an  attractive  method  for  experimentally 
determining  plasma  parameters  in  both  laboratory  and  space  plasmas.  Recently,  several  authors^*^ 
have  employed  classical  Langmuir^  (single)  probes  in  the  characterization  of  low  power  arcjet 
plumes.  While  single  probes  are  relatively  simple  to  construct,  it  is  necessary  to  sweep  the  bias 
voltage  supplied  to  the  probe  while  measuring  the  V-I  characteristic.  In  addition,  it  is  a  significant 
effort  to  extract  the  electron  density  11*  and  electron  temperature  Tc  from  the  corresponding  V-I 
characteristic.  Alternatively,  a  triple  probe  system^  can  be  employed  which  allows  sir^ultaneous 
measurement  of  Te  and  tie  without  the  necessity  of  a  voltage  sweep.  The  feasibility  of  applying 
triple  probes  to  electric  propulsion  devices  has  been  demonstrated  in  MPD  thruster  plumes^  and, 
more  recently,  in  low  power  arcjet  plumes.^  In  addition,  Habiger,  et  al.^  have  recently  presented  a 
general  review  of  electrostatic  probe  techniques  and  their  uses  in  electric  propulsion  devices. 

Using  the  theory  developed  originally  by  Kanal,*®  several  authors  have  measured  velocities 
in  flowing  plasmas  using  a  crossed-probe  technique. *2  By  combining  the  crossed-probe 
method  and  the  triple  probe  technique  discussed  above,  a  quadruple  electrostatic  probe  has  been 
devised  and  its  use  demonstrated  in  an  MPD  thruster. The  quadruple  probe  is  capable  of 
simultaneously  yielding  T^,  ne,  and  the  ratio  of  the  plasma  velocity  to  the  most  probable  thermal 
speed  (Uj/Cm).  The  goal  of  the  present  research  is  to  demonstrate  the  usefulness  of  the  quadruple 
probe  in  the  near-field  plume,  less  than  1.5  mm  from  the  thruster  exit  plane,  of  a  low  power  arcjet 
thruster. 

Quadruple  Probe  Schematic  and  Potential  Plot 

Figures  2.1a  and  2.1b  show  an  electrical  schematic  and  potential  plot  of  the  quadruple 
probe,  consisting  of  three  cylindrical  probes  P],  P2,  P3  which  are  aligned  with  the  plasma  flow 
vector  and  one  probe  P4  which  is  perpendicular  to  the  plasma  flow.  P3  and  P4  are  biased  at 
constant  voltages  Vd3  and  V(i4  relative  to  Pi,  and  P2  assumes  the  floating  potential  Vf  of  the  local 
plasma  (l2=0).  Probes  1,  3  and  4  are  also  electrically  floating  so  that  ion  current  collected  by  P3 
and  P4  is  balanced  by  the  electron  current  collected  by  Pj .  Three  quantities  are  recorded:  currents 
I3  and  I4  and  the  potential  difference  V(j2. 


*Ph.D.  Candidate,  Scott  A.  Bufton,  leads  our  group  in  advanced  probing  in  the  exit  plane  of  a  1  kW  class  arcjet. 
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Fig.  2.1a.  Quadniple  probe  electrical  schematic 
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Fig.  2.1b.  Quadruple  probe  potential  plot  for  constant  plasma  potential  Vp. 
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The  addition  of  the  perpendicular  probe  P4  allows  the  determination  of  the  parameter  Uj/Cm, 
determined  from  the  ratio  I4/I3  for  crossed  probes  with  =  Vd4.  Currents  I3  and  I4  are 
determined  by  measuring  voltage  drops  across  resistois.  In  Fig.  2.1,  voltages  Vd3  and  Vd4  are 
provided  by  batteries,  with  Vd3=Vd4=12V. 

Quadruple  Probe  Theory 

Quadruple  probe  analysis  assumes:  1)  The  ion  sheath  surrounding  the  tungsten  electrodes 
is  thin  compared  to  the  wire  radius  rp-,  2)  Each  probe  wire  and  sheath  are  collisionless  (>^fp  »  rp 

»  sheath  thickness);  3)  The  electron  energy  distribution  is  Maxwellian;  and  4)  A  wake  region 
of  low  ion  and  electron  density  downstream  of  P4  generates  collection  areas  for  ions  and  electrons 
of  A4i=A/7i  and  K^=kJl,  respectively. 

Derivation  of  quadruple  probe  response  is  similar  to  that  of  the  triple  probe^  and  has  been 
presented  elsewhere. the  present  study,  the  probe  is  applied  to  a  simulated  hydrazine 
plume,  with  multiple  species  of  ions  present,  whereas  previous  applications  of  triple  and  quadruple 
probes  have  been  for  plasmas  with  one  dominant  species  of  positive  ion.  Single  Langmuir  probes 
have  been  used  in  nitrogen  and  hydrogen-based  arcjet  plumes,'-^  but  their  use  is  not  contingent  on 
knowledge  of  the  species. 

For  a  quadruple  probe  with  Ai=A2=A3=A,  A4i=A/n,  A4e=A/2  and  Vd3=Vd4,  the  following 
expression  can  be  derived: 

1  - 1  +  exp(<|>Vd3)  -  2  exp  [<|>  (Vd3~Vd2)] 

jlexp[(t>(Vd3-Vd2)]-i)  (2.1) 

Since  Vd3  is  prescribed  and  Vd2  is  experimentally  determined,  Eq.  2.1  determines  Tg  =e/k<!>  through 
iteration.  Knowledge  of  the  plume  species  is  not  required  to  determine  electron  temperature. 

The  ion  saturation  current  density  jj  is  related  to  measured  quantities  by 

ji= _ Ui _ 

2A[exp(<|>Vd2)-(^+i)]  (2-^) 

Normally,  Eq.  2.2  is  used  to  determine  ne  from  I3  and  Vd2  by  relating  jj  to  the  electron  density 
through  the  Rohm  sheath  analysis.  For  a  plasma  with  only  one  species  of  positive  ion,  the  Rohm 
analysis  gives 
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(2.3) 


For  a  plasma  with  multiple  ion  species  Eq.  2.3  must  be  modified.  Assuming  that  each  collected 
positive  ion  behaves  independently  as  it  accelerates  through  the  Bohm  sheath  toward  the  probe,  jj 
becomes: 


'#  species 
I 

(-47-)] 

2 

.  j=l 

l(inj)^/2| 

(2.4) 


Assuming  equilibrium  dissociation  of  2H2  +  N2  in  the  arcjet  constrictor  with  each  species 
concentration  frozen  between  the  constrictor  and  the  exit  plane,  numerical  equilibrium  codes  predict 
H+  and  N+  as  the  only  positive  ions  present  in  measurable  levels, *5  with  0.4  <  nN+ZuH"^  <  0.7  for 
most  equilibrium  temperatures  of  interest  (Fig.  2.2).  The  Bohm  ion  saturation  current  is  then: 


ji  =  e(kTe)l/2exp(-l) 


'  npp-  ,  nM»  \ 
(mH+)^''^  (mN+)^''^/ 


(2.5) 


For  n€=nH'*'  +  nN"*",  Eqs.  2.2  and  2.5  yield 

K(3l3exp(l/2)(mH-^)^/^) 

2  A  e  (kTe)l/2  Jexp(^Vd2)  *  (UJ-)] 


(2.6) 


where  K  is  the  correction  to  the  electron  density  due  to  multiple  ion  species: 


1+M 
_ W* 

1  +  1 

jl/2 

(2.7) 


Note  that  Eq.  2.6  neglects  velocity  slip  between  each  species  of  ion  and  radial  diffusion  of 
H+  and  N+.  For  a  plume  containing  only  H+,  nN+  =  0.0,  k  reduces  to  unity,  and  Eq.  2.6  reduces 
to  single  component  form.  Figure  2.3  gives  k  as  a  function  of  the  ion  number  density  ratio  for  an 
H-N  plasma.  For  6,000  <  T  <  20,000  K,  Figs.  2.2  and  2.3  predict  1  <  K  <  1 .4. 
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Number  Density  Ratio, 


Fig.  2.3. 


Correction  factor  for  electron  density  ne  in  an  H-N  plasma  due  to  N"*’  ions. 


The  use  of  crossed  electrostatic  probes  for  measuring  plasma  velocity’ 0.12  the 
implementation  of  the  crossed  probe  technique  into  the  quadruple  probe’^  ’'*  have  been  presented 
elsewhere.  For  a  cylindrical  electrode  oriented  at  an  angle  0  relative  to  the  plasma  flow,  the 
collected  ion  current  is  a  function  of  the  quantity  (Ui/cni)sin0.  For  the  case  of  two  equally  biased 
probes,  with  one  probe  aligned  with  the  flow  vector  and  the  other  one  normal  to  it,  the  single 
species  collected  ion  current  ratio  is:’^ 


I//  I3  vs  A3 


exp(ui/Cm)^  I 

n  =  0 


(uj/  Cm)'' 


n! 


'Ml 


(2.8) 


Because  of  the  wake  effect  described  above,  the  current  collection  area  ratio  for  ions  is  A4/A3  = 
l/7t.  Combining  Eq.  2.8  with  measurements  of  the  currents  I3  and  I4  yields  the  quantity  Uj/Cm, 
which  is  related  to  the  plasma  Mach  number  by  M=  (2/Y)’^2(Ui/Cm),  where  y  is  the  specific  heat 
ratio  of  the  plasma.  If  information  about  the  ion  temperature  Tj  is  known  (Cm=(2kTi/mi)’^2)  the 
plasma  streaming  velocity  can  be  determined  from  Eq.  2.8. 

Apparatus  /  Experimental  Approach 

Electrostatic  probe  experiments  were  carried  out  at  the  exit  plane  of  a  1  kW  NASA  arcjet 
thruster^^  operating  at  a  flow  rate  of  50  mg/s  of  2H2  +  N2  to  simulate  fully  decomposed  hydrazine. 
Each  component  of  the  propellant  was  individually  metered  by  Unit  mass  flow  controllers  and 
mixed  in  the  propellant  line.  The  thruster  operated  at  10  A  and  104  to  1 10  volts,  depending  on 
how  long  the  arcjet  had  been  operating,  and  exhausted  into  a  1  m  diameter  x  1.5  m  long  vacuum 
tank  pumped  by  a  2500  cfm  Roots-type  blower.  Backing  for  this  was  provided  by  a  second  1300 
cfm  Roots  blower  and  two  150  cfm  displacement  pumps.  The  tank  pressure  during  arcjet 
operation  was  -200  mTorr. 

A  schematic  of  the  quadruple  probe  utilized  in  this  study  is  shown  in  Fig.  2.4.  The  probe 
is  comprised  of  4  individual  tungsten  wire  tips  of  0.25  mm  diameter  with  an  exposed  length  of  3.5 
mm.  Each  wire  mounts  in  round  four-bore  *dumina  (AI2O3)  tubing,  providing  rigidity,  electrical 
insulation,  and  high  heat  capacity  to  allow  probe  measurements  at  the  exit  plane  of  the  thruster. 
The  alumina  is  2.4  mm  in  diameter,  with  four  0.5  mm  diameter  through  holes  (Fig.  2.4),  and  is 
mounted  in  thin-walled  stainless  steel  tubing  to  provide  noise  shielding  and  to  facilitate  mounting 
the  probe  and  its  leads  to  the  positioning  mechanism.  In  early  experiments,  epoxy  was  used  to 
secure  each  wire  in  the  alumina,  but  the  epoxy  tended  to  ablate  away  over  time,  exposing  more  of 
each  tungsten  electrode.  This  method  was  abandoned  in  favor  of  a  high  temperature  ceramic-based 
adhesive  to  bond  each  wire  to  the  alumina. 
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Fig.  2.4.  Schematic  of  quadruple  probe  tip  swept  through  the  arcjet  plume. 
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The  thin-wall  stainless  steel  tube  is  rigidly  connected  to  two  ceramic  thermocouple  plugs 
with  a  compression  fitting  (Fig.  2.5).  Each  plug  houses  connections  for  two  tungsten  wires, 
allowing  connection  to  the  probe  electronics  by  a  single  four-prong  plug  located  at  the  rear  and  a 
vacuum  feedthrough  in  the  tank  wall. 

The  complete  assembly  is  mounted  to  a  stepper  motor  mounted  on  a  linear  translation 
carriage  (Fig.  2.5).  The  stepper  motor  allows  the  probe  axis  to  be  rotated  in  increments  of  0.9°. 
The  probe  arm  is  coupled  directly  to  the  motor  shaft  with  the  probe  sensing  tip  on  the  shaft  axis, 
allowing  probe  angle  to  be  varied  without  changing  the  tip  location.  For  the  present  study,  the 
probe  axis  is  always  parallel  to  the  thruster  axis. 

The  motor  and  probe  assembly  is  mounted  to  a  carriage  capable  of  a  +  5  cm  linear 
translation  perpendicular  to  the  arcjet  axis,  allowing  the  quadruple  probe  assembly  to  be  manually 
swept  through  the  plume.  The  carriage  platform  is  mounted  to  a  precision  shaft/linear  ball  bearing 
assembly  giving  smooth  and  repeatable  probe  movement  by  a  shaft  passing  through  a  vacuum 
feedthrough  in  the  tank  wall.  Carriage  and  probe  position  are  determined  by  a  precision  10  k£2 
linear  potentiometer  coupled  to  the  carriage  shaft.  The  probe  position  is  determined  to  0.05  mm 
relative  to  the  arcjet  nozzle  by  calibrating  the  10  kQ  linear  potentiometer  output  against  the  outer 
diameter  of  the  nozzle.  The  probe  system  can  be  removed  as  a  unit  from  the  vacuum  system  and 
bench-tested. 

Quadruple  Probe  Preliminary  Results  and  Analysis 

Preliminary  experimental  results  are  presented  for  the  near-fteld  plume  of  the  thruster,  with 
the  probe  tip  located  1.5  mm  from  the  thruster  exit.  Electron  densities  are  calculated  using  Eq.  2.6, 
with  a  correction  factor  k=1.35,  which  corresponds  to  nN'^/nH'^=0.55.  Variation  of 
0.4<nN'^/nH'^<0.7  causes  only  a  ±  5%  error  in  k  and,  hence,  in  tie.  Electron  number  density  results 
are  shown  in  Fig.  2.6  for  probe  sweeps  recorded  at  -4  minute  intervals  during  arcjet  warm-up. 
Low  Me  is  measured  by  the  probe  outside  of  the  exit  plane  diameter,  as  expected.  The  electron 
temperature  used  to  determine  nc  is  1  eV  for  each  case,  as  discussed  below.  For  the  range  0.3  <  Te 
<  1.0  expected  at  the  exit  plane,  ne  varies  approximately  as  Te"!^^,  so  that  a  corresponding 
correction  to  the  ne  values  of  Fig.  2.6  must  be  made.  The  increases  in  ne  with  time  in  Fig.  2.6  are 
not  due  to  variations  in  Te.  since  a  factor  of  3  decrease  in  Te  would  be  required  to  explain  the 
measured  increase  in  ne. 

The  thruster  appears  to  be  operating  in  a  steady  state  condition,  based  on  probe 
measurements,  after  approximately  4  to  8  minutes.  This  is  consistent  with  previously  reported 
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Fig.  2.5.  Schematic  of  quadruple  probe  mount.  Probe  axis  is  parallel  to  the  thmster 
in  this  study. 


Radius,  mm 

Fig.  2.6.  Electron  density  ng  versus  radius  during-thruster  warm-up. 
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data  indicating  that  the  arcjet  operating  voltage  level  is  98%  of  its  steady  state  level  after 
approximately  5  minutes  of  operation. 

Since  the  quadruple  probe  is  aligned  with  the  thruster  axis  during  these  sweeps,  a 
correction  is  required  to  account  for  the  effect  of  flow  transverse  to  the  probe  axis.  For  probe 
locations  near  the  plume  centerline,  these  corrections  are  small.  At  locations  far  from  the  thruster 
axis  this  correction  may  be  quite  large,  but  the  thin  sheath  assumption  breaks  down  in  these 
regions  and  a  correction  was  not  attempted.  Assuming  a  correction  to  the  !)«  values  in  Fig.  2.6 
consistent  with  a  centerline  Te  of  0.5  eV,  our  data  are  lower  than  the  ne  results  of  other 
researchers*"^  by  a  factor  of  ~2  -  4. 

The  measured  currents  I3  and  I4  and  the  current  ratio  I4/I3  used  to  determine  Ui/Cm  from  Eq. 
2.8  are  plotted  in  Fig.  2.7  for  t~8  minutes  after  arcjet  start-up.  For  r<3  mm  near  the  centerline, 
I4/I3  -  0.6,  corresponding  to  Ui/Cni=l.S  and  a  Mach  number  M~1.9.  This  calculation  neglects  the 
effect  of  multiple  species  on  the  ion  current  to  the  crossed  probes.  For  dissociated  hydrazine  at  0.5 
eV,  with  a  mean  molecular  weight  of  5.3,  the  sound  speed  is  a  -  3,300  m/s,  and  the  corresponding 
M=1.9  flow  velocity  is  6,200  m/s.  This  low  Mach  number  is  consistent  with  results  for  high  area 
ratio,  low  Reynolds  number  nozzles. 

A  typical  V(i2  trace  versus  position  across  the  thruster  face  is  shown  in  Fig.  2.8.  Since,  by 
Eq.  2.1,  Te  increases  monotonically  with  Vd2,  the  non-symmetric  trace  would  seem  to  imply 
the  unexpected  result  of  non-symmetric  Te,  with  Te  increasing  monotonically  from  a  minimum  of 
-0.5  eV  at  r=3  mm  to  -0.7  eV  on  axis  to  -1.2  eV  at  the  right  side  of  the  exit  plane.  Because  of  the 
symmetry  in  the  De  profiles,  this  amount  of  Te  asynunetry  seems  unlikely. 

Several  experiments  were  performed  to  investigate  the  nature  of  the  observed  Vd2 
asymmetry.  Probe  heating  and  thermionic  emission  of  electrons  were  first  eliminated  by 
examining  Vd2  traces  obtained  by  passing  the  probe  through  the  plume  in  each  of  two  possible 
directions.  The  resulting  traces  in  this  case  were  mirror  images  of  each  other,  indicating  that  probe 
heating  was  not  a  factor.  To  rule  out  physical  asymmetry  in  the  plume,  the  entire  thruster  was 
rotated  180*  about  the  thrust  axis.  This  experiment  yielded  results  identical  to  those  recorded 
before  the  arcjet  was  rotated. 

Probe  asymmetry  itself  was  then  investigated.  Initial  experiments  were  conducted  with  a 
quadruple  probe  configured  as  shown  in  Fig.  2.9a.  A  more  symmetric  configuration.  Fig.  2.9b, 
showed  improvement  but  not  elimination  of  the  asymmetry. 
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Fig.  2.7.  Probe  currents  I3, 14,  and  I4/I3  giving  M~  1 .9  in  the  center  of  the  plume. 


Fig.  2.8.  Typical  Vd2  vs.  r  showing  asymmetry  due  to  Vp  gradients,  and  Te~0.6  eV 
on  axis. 
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Current  Ratio,  1  / 1 


Fig.  2.9.  Quadruple  probe  configurations  used  to  change  probe  syrmnetry  (a,  b)  and 
for  radial  E-field  measurements  (c). 
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It  thus  appears  likely  that  radial  gradients  in  the  plasma  parameters  (Te.  ne,  Vp)  across  the 
probe  volume  cause  the  asymmetry.  To  verify  this,  the  double  floating  probe  configuration  shown 
in  Fig.  2.9c  was  used  to  measure  the  radial  electric  field  Ef.  For  two  floating  probes  in  a  plasma 
with  sheath  potential  Vj  and  a  potential  difference  between  the  two  electrodes  AVf  =V(7  -Vfi,  the 
difference  in  plasma  potential  is: 


Vp2-  Vpi  =  Vf2-  Vfi  - (Vs2-  Vsi)  (2.9) 

The  results  of  a  AVf  measurement  are  plotted  in  Fig.  2. 10.  The  difference  in  sign  between 
the  two  peaks  at  r=S.S  mm  is  a  consequence  of  electrode  orientation  remaining  constant  while  the 
plasma  gradients  change  sign  from  one  side  of  the  plasma  plume  to  the  other.  From  Eq.  2.9,  the 
change  in  plasma  potential  across  the  probe  can  be  found  if  AVf  and  the  probe  sheath  drops  are 
known.  For  a  single  floating  electrode,  the  sheath  potential  drop  is  given  in  the  usual  way  by 
summing  the  ion  and  electron  current  into  the  probe  and  solving  for  the  potential  in  the  Boltznuum 
term  that  satisfies  the  zero  net  current  requirement.  This  yields  a  sheath  potential  of 


Assuming  that  the  logarithm  term  does  not  vary  much  over  the  probe  separation, 

Vs2-Vsi=AVs  =  ;^ATeln(“Ly  (2.11) 

2e  \“'eTi/ 

If  we  assume  that  the  gradients  in  electron  temperature  are  small  over  the  electrode  separation,  then 
AVs=0  and  measurements  of  AVf  will  equal  the  plasma  potential  difference  Vp2-Vpi,  and  radial 
electric  field  Er  =AVf(d  at  the  thruster  exit.  The  resulting  magnitude  of  Ef  versus  position  is  shown 
in  Fig.  2.1 1.  The  field  is  approximately  zero  at  the  center  as  expected,  with  peaks  of  ~  200  V/m 
near  the  edges  of  the  thruster  exit  plane.  From  Er  =  -dVp/dr,  we  can  integrate  the  electric  field  to 
give  the  profile  of  Vp  across  the  thruster  exit  (Fig.  2.12)  showing  a  plasma  potential  profile  with  a 
shallow  potential  well  of  depth  -0.1  V,  with  its  minimum  near  the  thruster  axis  and  steep  gradients 
outside  the  exit  plane  at  r  >  5  mm. 

The  probe  measurements  are  strong  evidence  for  the  existence  of  current  density  and  ohmic 
heating  in  the  near-field  plume.  From  the  mass  flow  rate  and  measured  tie,  the  estimated  ionization 
fraction  on  the  centerline  is  a*  10-^,  giving  a  plasma  electrical  conductivity  of  a«100  mho/m.  In 
the  exit  plane  Er  -  100  V/m,  equivalent  to  j=aE~l  A/cm^  crossing  the  exit  plane  of 
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Fig.  2.10.  V2f-Vif  measurement  with  E-field  configuration. 


Radius,  mm 


Fig.  2.11.  Radial  E-field  versus  position  across  thruster  face. 
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the  thruster  and  reattaching  at  the  comer  of  the  nozzle.  If  this  current  is  flowing  into  a  region  of 
approximately  0.2  cm^,  the  total  current  crossing  the  exit  plane  amounts  to  ~0.2  A,  or  2%  of  the 
total  current  supplied  to  the  arcjet.  Furthermore,  this  current  represents  a  small  amount  of  ohmic 
heating  (j^/<J~l  W/cm^)  in  the  exit  plane  region. 

The  indication  seen  here  of  current  in  the  near-field  plume  is  consistent  with  the  condition 
of  j  perpendicular  to  the  nozzle  wall,  which  would  force  current  into  the  plume.  These  preliminary 
measurements  also  may  carry  implications  for  choosing  downstream  boundary  conditions  for  arcjet 
numerical  models. 

Implications  of  Non-uniform  Plasma  Potential 

The  basic  triple  and  quadruple  probe  theories^>i^  assume  zero  AVp,  strictly  tme  here  only 
on  axis.  Figure  2.13  is  a  potential  plot  for  the  quadruple  probe  in  a  plasma  potential  gradient  (E- 
field).  Because  biasing  voltages  and  Vd4  are  large,  small  changes  in  Vp  have  little  effect  on  the 
ion  saturation  current.  However,  the  gradient  in  Vp  can  affect  the  measured  Te-  It  can  be  shown 
that  the  quadruple  probe  (qp)  voltage  Vd2  (Fig.  2.13)  is  given  more  generally  by: 


Vd2  =  IVim-Vlmlqp  -  (Vp2  -  Vpi)  (2.12) 

Substituting  for  Vp2-Vpi,  as  measured  by  a  double  floating  probe  (dfp),  from  Eq.  2.9: 

Vd2=IV2ni  -V  Imlqp  -[(V  f2*Vf|)dfp  -(V  s2-Vsl)dfp]  (2. 1 3) 

Combining  Eqs.  2.1 1  and  2.13,  the  electron  temperature  can,  in  principle,  be  determined  despite 
the  existence  of  gradients  in  Vp.  This  work  is  in  progress. 
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Fig.  2. 


Potential  plot  for  the  quadruple  probe  in  a  plasma  potential  gradient  (E-field), 
where  Vp];6Vp2. 
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2.2  Quadruple  Probe  Modeling  with  Plasma  Gradients 
Introduction 

Preliminary  results  of  quadruple  probe  measurements  at  the  thruster  exit  plane  show 
smooth  symmetric  profiles  for  ne  (Fig.  2.6)  and  the  electrode  currents  I3  and  I4  (Fig.  2.7). 
However,  Fig.  2.8  shows  an  asymmetric  V(i2  profile,  which  directly  implies  an  unexpected 
asynunetry  in  the  measured  Te  profile  through  Eq.  2.1.  Several  possible  explanations  for  this 
asymmetry  have  been  considered,  including:  1)  electrode  heating  and  thermionic  emission,  2)  the 
presence  of  a  physical  asymmetry  in  the  plume,  3)  probe  geometry,  and  4)  the  influence  of 
gradients  in  plasma  parameters  over  the  probe  volume.  Each  of  the  first  three  reasons  were 
eliminated  as  possibilities  through  a  judicious  choice  of  experiments,  leaving  the  fourth  explanation 
as  a  remaining  possibility. 

The  quadruple  probe  theory  derived  in  Section  2.1  assumes  that  the  various  plasma 
parameters  (Te,  ne  and  Vp)  are  uniform  over  the  probe  volume.  It  has  been  hypothesized,  despite 
the  small  (~  1  mm)  width  of  the  quadruple  probe,  that  slight  gradients  in  the  plasma  parameters 
over  the  probe  volume  may  be  responsible  for  the  asymmetric  Te  proEles  measured  at  the  exit 
plane.  This  possibility  is  being  investigated  by  examining  the  effect  of  selected  radial  profiles  of  Te 
and  rie  on  the  quadruple  probe  results  when  the  current  to  each  electrode  is  evaluated  at  the  local  Te 
and  lie  of  each  electrode.  More  precisely,  a  quadruple  probe  computational  model  has  been 
constructed  that  effectively  runs  the  quadruple  probe  experiment  in  reverse.  In  this  way,  one  can 
specify  profiles  of  Te  and  tie  and  calculate  the  quantities  normally  measured  by  the  quadruple 
probe:  the  difference  in  the  potentials  of  electrodes  1  and  2  (V(]2);  and  the  measured  current  drawn 
by  electrode  3  (I3).  The  predicted  values  of  V(i2  and  I3  can  then  be  compared  to  the  corresponding 
measured  values  to  gain  insight  into  the  extent  that  the  finite  dimensions  of  the  electrostatic  probe 
influence  the  measured  profiles  of  Te  and  ne.  Ultimately,  the  model  may  provide  a  means  of 
solving  for  the  Te  and  tie  profiles  from  the  measured  asymmetry  in  the  Vd2  signal. 

Model  Derivation 

In  the  laboratory,  the  measured  quadruple  probe  parameters  Vd2  I3  are  related  to  Te  and 
ne  through  Eqs.  2.1  and  2.6.  To  develop  the  quadruple  probe  model  in  the  presence  of  plasma 
gadients,  the  four  individual  electrode  current  equations  are  rewritten  in  terms  of  the  plasma 
parameters  near  each  electrode,  instead  of  assuming  uniformity  in  Te  and  ite  across  the  probe 
volume.  The  electrode  current  equations  are  rearranged  to  solve  explicitly  for  Vd2  and  I3  in  terms 
of  specified  profiles  of  Te  and  He-  This  results  in  the  following  relations  for  the  performance  of  the 
quadruple  probe: 
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V2  =  V2  (Te2.  ne2) 

(2.14) 

Vi  =  Vi  (Tei,Te3,Te4.nel.  ne3.  ne4) 

(2.15) 

Vd2  =  V2-V, 

(2.16) 

I3  =  I3  (Te3.  ne3) 

(2.17) 

The  subscripts  on  the  parameters  Te  and  ne  in  the  above  equations  denote  the  local  value  of  those 
parameters  at  each  of  the  four  electrodes  (1-4)  making  up  the  quadruple  probe. 

Initial  versions  of  this  model  have  been  implemented  in  a  spreadsheet  format  that  solves  for 
V(i2(r)  and  l3(r),  for  given  profiles  of  Te(r)  and  ne(r).  Since  the  centerline  experimental  data  are,  by 
symmetry,  at  zero  gradient  conditions,  these  data  are  used  to  determine  the  centerline  values  of  Te 
and  ne.  For  the  initial  version  of  the  model,  Te(r)  and  ne(r)  are  specified  as  Gaussian  distributions, 
with  a  peak  value  corresponding  to  the  experimentally  measured  thruster  centerline  values.  The 
width  (standard  deviation)  of  each  profile  is  adjusted  to  study  the  effect  on  the  predicted  quantities 
V(i2  and  I3.  '^ince  (Eq.  2.16)  and  I3  (Eq.  2.17)  are  strong  functions  of  Te  and  ne,  respectively, 
there  is  sufficient  uncoupling  that  the  model  can  be  used  to  solve  iteratively  for  the  Te(r)  and  ne(r) 
profiles  to  reproduce  the  measured  Vd2  and  I3  profiles.  The  numerical  program  moves  the  probe 
position  across  the  thruster  face  in  -0.1  mm  increments  (-1/10  of  the  probe  width)  while  solving 
for  the  predicted  values  Vd2  and  I3,  which  are  then  compared  to  the  measured  values  Vd2  and  I3. 
The  Te  and  ne  profiles  are  then  adjusted  to  match  the  measured  values  as  the  process  is  repeated. 

Probe  Model  Preliminary  Results 

Figures  2.14-2.16  show  preliminary  results  of  the  above  analysis  for  one  set  of 
experimental  data.  Figure  2.14  shows  the  Gaussian  profiles  of  Te(r)  and  ne(r)  that  best  reproduced 
the  measured  profiles  of  V(i2  and  I3,  along  with  the  physical  boundary  of  the  exit  plane  and  the 
thruster  centerline.  The  width  of  ne(r)  corresponds  to  a  standard  deviation  (o)  of  3.0  mm,  while 
the  Te(r)  profile  is  much  broader  with  a=8.0  mm.  Figures  2.15  and  2.16  compare  the  measured 
and  resulting  predicted  values  of  I3  and  Vd2,  respectively.  These  figures  show  good  agreement  in 
trends  between  the  measured  and  predicted  values,  which  is  strong  evidence  that  the  quadruple 
probe  measurements  are  being  influenced  by  gradients  in  plasma  parameters  over  the  probe 
volume.  The  agreement  is  poorest  far  from  the  centerline,  presumably  because  basic  assumptions 
in  the  quadruple  probe  theory  break  down  in  this  region.  In  particular,  the 
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Fig.  2. 14.  Gaussian  radial  profiles  of  Te  and  ne  at  the  thruster  exit  plane  extracted  from  the 
measured  values  of  I3  and  Vd2. 
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Fig.  2.  IS.  Comparison  of  the  measured  I3  to  the  value  predicted  by  the  quadruple  probe  model 
for  Te(r)  and  ne(r)  with  a  of  8.0  and  3.0,  respectively. 
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Fig.  2.16.  Comparison  of  the  measured  Vd2  to  the  value  predicted  by  the  quadruple  probe 
model  for  Te(r)  and  nc(r)  with  ct  of  8.0  and  3.0,  respectively. 
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parallel  probes  are  misaligned  with  the  flow  vector  and  the  Debye  sheaths  around  each  electrode 
become  thick  and  interact  with  each  other  far  from  the  thruster  centerline. 

Conclusions 

Electrode  sheath  constraints  prevent  the  construction  of  electrostatic  probes  that  are  small 
enough  to  avoid  the  problem  of  gradients  of  Te  and  ne  over  the  probe  size.  In  order  to  continue 
using  this  low-cost  and  robust  measurement  technique  in  small  (large  gradient)  arcjet  devices,  it  is 
necessary  to  understand  how  these  gradients  affect  probe  performance.  In  addition,  it  is  necessary 
and  desirable  to  be  able  to  extract  the  profiles  of  Te  and  tie  from  the  raw  data  (Fig.  2.8).  To  this 
end,  the  gradient  model  appears  capable  of  replicating  the  asynunetric  measured  probe  profiles 
inside  the  nozzle  exit  (0<r<5  mm).  Future  versions  of  the  model  are  being  developed  which  will 
optimize  the  agreement  between  the  measured  and  predicted  values  by  varying  the  shape  of  the  Te 
and  ne  profiles. 
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2.3  Electrostatic  Probe  Velocimetry 
Overview 

In  addition  to  the  quantities  ne  and  Te,  quadruple  probe  measurements  also  yield  the 
quantity  uj/Cn,  (=  ion  velocity  /  most  probable  thermal  speed).  Since  Cni=(2kTi/mi)J''2,  knowledge 
of  the  heavy  particle  velocity  Ui  allows  determination  of  the  heavy  particle  temperature  Ti. 
Recently,  a  low-cost  method  of  measuring  velocities  in  arcjet  plumes  has  been  developed  and 
demonstrated  by  Spores,  et  al.  on  1  kW  hydrogen  arcjets.®'^*  In  this  method,  called  Current 
Modulation  Velocimetry  (CMV),  a  short  duration  current  spike  is  superimposed  on  the  arcjet 
operating  current.  The  effect  of  this  current  spike  on  the  plasma  is  then  monitored  in  a  time-of- 
flight  fashion  as  it  is  convected  past  two  measurement  positions  in  the  thruster  plume.  This 
research  will  utilize  a  modified  and  improved  version  of  the  CMV  technique,  coupled  with 
quadruple  electrostatic  probe  measurements,  to  attempt  to  measure  radial  profiles  of  Te,  ne,  Uj  and 
Tj  in  the  near-field  plume  of  a  1  kW  arcjet. 

Review  of  the  CMV  Technique 

In  the  work  by  Spores,  et  al,  a  short  duration  current  pulse  was  superimposed  on  the  arcjet 
operating  current.**  1*  Figure  2.17  shows  the  effect  of  one  such  current  pulse  on  the  arcjet  current 
and  voltage.  Also  shown  in  Fig.  2. 17  is  the  periodic  ripple  associated  with  the  arcjet  current  and 
voltage,  which  is  a  consequence  of  the  switching  nature  of  the  arcjet  power  processing  unit  (PPU). 
The  current  and  voltage  ripple  are  180‘’  out  of  phase,  which  is  a  result  of  the  negative  impedance 
characteristic  of  the  arcjet.  The  current  pulse  is  initiated  by  momentarily  closing  a  switch  which 
introduces  an  RC  circuit  parallel  to  the  arcjet.  Upon  switch  closure,  the  arcjet  voltage  appears 
across  the  RC  circuit,  which  charges  over  a  few  1-ps  time  constants.  Since  the  current  from  the 
PPU  is  held  constant  over  this  time  interval,  the  charging  capacitor  represents  a  short-lived  current 
deficit  in  the  arcjet.  The  arcjet  current  then  rings  sharply  upward,  significantly  passing  its  previous 
maximum  before  it  quickly  damps  to  its  normal  PPU  ripple. 

It  is  believed  that  the  positive  ring  of  the  current  spike  causes  a  momentary  increase  in 
electron  recombination,  which  repopulates  the  uppermost  energy  levels.  These  electrons  cascade 
downward,  giving  off  radiation  consistent  with  each  electron  transition.  The  propellant  in  the  arc¬ 
heating  region  during  the  pulse  is  therefore  "tagged."  Velocity  is  inferred  by  measuring  the  time 
required  for  this  tagged  flow  to  traverse  the  distance  between  two  observation  points.  In  Refs.  8  & 
18,  this  effect  was  monitored  by  measuring  changes  in  H  Balmer  radiation  at  two  different  axial 
locations  using  photomultiplier  tubes  (PMT),  Hydrogen  Balmer-alpha  notch  filters 
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Fig.  2. 17.  Arcjet  current  and  voltage  traces  showing  the  current  pulse  used  in  Refs.  8  &18 

and  the  periodic  ripple  associated  with  the  arcjet  PPU. 

(Taken  from  Spores,  et  al.^*) 
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and  the  appropriate  optics.  The  delay  between  the  emission  signals  measured  by  each  PMT  is 
determined  and  related  to  velocity  through  the  known  detector  control  volume  separation. 

While  optically  monitoring  the  effect  of  the  electron  recombination  usually  proved  useful  in 
measuring  plasma  velocities,  it  was  noted  that  the  positive  ring  of  the  current  modulation  d  not 
always  trigger  the  electron  cascade  phenomenon  necessary  to  cause  the  requisite  Balmer  emission. 
In  fact,  it  was  found  that  for  hydrogen  arcjet  operating  voltages  below  a  certain  minimum  voltage, 
the  CMV  technique  described  above  did  not  work. 

Modifications  and  Improvements  in  CMV  for  use  at  UIUC 

The  optical  (PMT)  measurement  technique  used  in  Refs.  8  &  18,  while  non-intrusive,  is  a 
line-of-sight  technique  that  precludes  the  direct  measurement  of  spatially  resolved  velocities.  The 
measured  velocity  is  a  radially  integrated  value  over  the  PMT  line-of-sight.  This  may  explain  the 
shape  of  the  PMT  signal  profiles  measured  by  Spores,  et  al.  Both  signals  rise  rather  steeply 
initially  ( -1  ps)  and  fall  off  much  more  gradually  (-6  ps).  The  sharp  initial  rise  is  possibly  due  to 
emission  from  the  higher  velocity  centerline  plasma,  while  both  flow  divergence  and  slower  off 
centerline  plasma  may  contribute  to  the  gradually  decreasing  segment  of  the  PMT  signals.  Ideally, 
a  narrow  width  signal  is  desired  for  this  type  of  measurement,  and  the  wide  shape  of  the  signals 
from  optical  methods  may  introduce  additional  errors  during  calculation  of  the  signal  time  delay. 

As  an  improvement  to  the  optical  measurement  methods  described  above  and  in  Refs.  8  & 
18,  a  time-of-flight  electrostatic  probe  system  will  be  utilized  to  measure  the  tagged  plasma  flow 
velocity.  The  electrostatic  probe  system  affords  several  advantages  over  optical  methods;  1)  It  is 
very  low-cost  and  relatively  robust,  2)  Measurements  can  be  spatially  resolved,  allowing 
measurements  of  velocity  profiles,  3)  Measured  signal  widths  are  likely  to  be  narrow,  making 
signal  delays  simpler  to  distinguish,  4)  Momentary  changes  in  basic  plasma  parameters  (Te,  ne, 
etc.)  can  be  monitored  to  indicate  passage  of  the  tagged  flow,  and  5)  Successful  velocity 
measurements  are  not  contingent  upon  causing  an  optical  event,  a  problem  encountered  by  other 
researchers.®’!* 

Time-of-Flight  Velocity  Probe  Preliminary  Design 

Figure  2.18  is  a  schematic  showing  a  preliminary  concept  for  the  time-of-flight  sensing 
probe  design.  The  probe,  similar  in  construction  to  the  quadruple  probe  used  for  ne,  Te  and  Uj/cm 
measurements,  consists  of  two  tungsten  electrodes  in  line  with  the  flow  vector  and  electrically 
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Fig.  2.18  Side  view  schematic  of  the  time-of-flight  electrostatic  probe  used  for  near¬ 
field  plume  velocity  measurements. 
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insulated  as  shown.  Although  this  probe  is  physically  intrusive,  appropriate  length  scales  indicate 
that  the  probe  electrodes  (0.25  mm  dia.)  will  be  in  the  free  molecular  flow  regime  and  are  thus 
expected  to  have  little  effect  on  the  flow.  Previous  experiments  with  the  perpendicular  electrode  of 
the  quadruple  probe  have  shown  that  an  exposed  electrode  length  of  only  a  few  millimeters 
provides  sufficient  collection  area  and  signal  strength.  The  two  electrodes  may  be  offset  slightly 
(Fig.  2.18)  to  avoid  an  attenuation  of  the  signal  at  electrode  2  due  to  the  possibility  of  charge 
depletion  caused  by  current  drawn  at  electrode  1,  the  upstream  electrode.  The  axial  probe 
separation  is  expected  to  be  approximately  S  mm  or  less.  However,  the  10  MHz  sampling  speed 
of  the  existing  data  acquisition  equipment  may  dictate  a  slightly  larger  distance  in  order  to 
accurately  resolve  the  time  delay  between  the  measured  fluctuations  at  each  electrode.  To  facilitate 
increasing  the  accuracy  of  the  time  delay  measurement  and  decreasing  the  electrode  axial 
separation,  faster  data  acquisition  equipment  will  be  examined. 

Monitoring  natural  fluctuations  in  the  ion  saturation  current  has  been  suggested^’^^  as  a 
means  of  using  time-of-flight  electrostatic  probes  to  measure  velocities  in  perturbed  plasma  flows. 
This  method  is  implemented,  but  instead  will  be  used  with  artificially  supplied  current  fluctuations. 
The  two  electrodes  (Fig.  2.18)  are  biased  sufficiendy  negative  with  respect  to  the  plasma  potential 
Vp  that  they  each  measure  the  local  ion  saturation  current.  From  a  particle  flux  standpoint,  the 
saturation  current  is  roughly  proportional  to  the  local  charged  particle  density  (and  also  collection 
area,  particle  velocity  and  particle  charge).  In  previous  work,  Sankovic^  showed  that,  for  a  1  kW 
hydrazine  arcjet,  the  electron  density  in  the  plume  (~10  cm  from  nozzle  exit)  was  a  strong  function 
of  the  thruster  specific  power.  More  specifically,  it  was  demonstrated  with  steady  state 
measurements  that  a  20%  decrease  in  arc  current  caused  nearly  a  30%  decrease  in  measured 
electron  density.  With  quasineutrality,  this  corresponds  to  a  similar  change  in  ion  saturation 
current.  This  is  to  be  expected  physically  since  changing  the  arc  current  changes  the  number  of 
energetic  electrons  available  for  ionizing  collisions.  Therefore,  superimposing  a  current  pulse  on 
the  steady  arcjet  operating  current  is  expected  to  cause  a  measurable  fluctuation  in  the  measured  ion 
saturation  current  as  the  tagged  flow  passes  each  electrode. 

Flow  velocity  will  be  calculated  using  the  time  delay  between  the  ion  saturation  current 
signals  measured  at  the  two  electrodes  and  the  known  electrode  separation.  The  method  of 
extracting  the  delay  betweei.  the  two  signals  is  contingent  on  the  shape  of  the  signals.  Two  very 
sharply  spiked  signals  would  require  little  more  than  a  visual  inspection  to  determine  the  time  lag. 
More  realistically,  the  signals  will  have  significant  widths  associated  with  them  and  the  process 
will  be  more  involved.  In  this  case  a  numerical  solution  may  be  required  to  determine  the  time 


54 


delay,  such  as  digitally  performing  a  Fast  Fourier  Transform  (FFT)  autocorrelation  routine  on  the 
two  signals. 

This  method  is  attractive  because  the  probe  electronics  can  be  rearranged  in  a  very  short 
time  to  measure  the  effect  of  the  current  pulse  shape  on  other  plasma  parameters.  The  polarity  of 
the  probe  bias  can  be  switched  so  that  the  electrodes  draw  electron  saturation  current,  which  will 
produce  a  much  stronger  signal  due  to  the  larger  mobility  of  the  electrons.  In  addition,  the  effect  of 
the  pulse  on  Te  can  be  examined  by  measuring  the  floating  potential  Vf  (relative  to  Vp)  of  each 
electrode,  since  Vf  varies  approximately  with  Te-  Finally,  a  double  probe  configuration  can  be 
used,  in  which  one  electrode  draws  ion  current  while  the  other  draws  an  equal  amount  of  electron 
current.  In  this  configuration  a  single  channel  can  be  monitored,  which  will  show  a  fluctuation 
each  time  the  tagged  flow  passes  an  electrode.  Since  changing  configurations  is  simple,  several 
can  be  attempted  to  check  for  agreement  or  to  determine  which  is  most  feasible  in  terms  of  signal 
quality  and/or  data  reduction. 

Introducing  the  Current  Spike  to  the  Arcjet  Circuit 

The  arcjet  electronic  circuit  is  shown  in  Fig.  2.19.  It  consists  of  three  components;  a 
NASA  LeRC  arcjet  power  processing  unit  (PPU),  a  1  kW  arcjet,  and  a  current  modulation  leg 
parallel  to  the  thruster.  The  current  modulation  section  of  the  circuit  is  connected  through  a 
MOSFET  gate  switch.  ^  **21  The  MOSFET  switch  and  its  controller  circuit  are  driven  by  a  pulse 
generator,  so  that  the  switch  closes  for  the  duration  of  the  pulse  supplied  by  the  pulse  generator. 
Therefore,  the  duration  of  the  current  modulation  ctm  be  varied  by  changing  the  width  of  the  pulse 
supplied  by  pulse  generator.  The  MOSFET  controlled  electrical  switch  shown  in  Fig.  2. 19  may  be 
used  in  conjunction  with  a  variety  of  methods  to  create  the  current  modulation  to  the  arcjet. 

A  few  variations  on  the  current  modulation  result  from  changing  the  components  in  the 
dashed  box  (Fig.  2.19),  as  shown  in  Figs.  2.20a  and  2.20b.  Figure  2.20a  shows  the  circuit  used 
by  Spores,  et  al.^®  Closing  the  MOSFET  switch  causes  a  momentary  current  loss  to  the  arcjet, 
followed  by  the  positive  current  "ringing"  spike.  It  is  this  positive  spike  that  causes  the  optical 
event  measured  in  the  above  mentioned  CMV  papers.  With  electrostatic  probes,  however,  the 
probes  need  only  sense  the  changes  in  Te,  ne  or  electrode  floating  voltage  Vf  as  the  tagged  flow 
passes  the  probes.  Therefore,  a  simple  resistive  load  may  be  sufficient  for  the  purposes  of  the 
present  research  (Fig.  2.20b).  A  momentary  closure  of  the  MOSFET  switch  will  cause  a  sudden 
and  short-lived  current  deficit  to  the  arcjet.  As  discussed  above,  this  should  cause  significant 
density  changes  such  that  the  effect  of  the  current  modulation  will  be  measurable  at  the  time-of- 
flight  electrostatic  probes. 
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Fig.  2.19  Electrical  circuit  schematic  for  MOSFBT  gate  controlled  modifications  of 
the  arcjet  current. 


Fig.  2.20  Electrical  circuit  elements  for  causing  the  arcjet  current  pulse/deficit. 
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Section  3:  Arcjet  Energy  Deposition  in  the  Constrictor  and  Nozzle 


3 . 1  Introduction 

Improvements  in  the  efficiency  and  performance  of  arcjet  thrusters  are  achieved  by 
understanding  the  fundamental  energy  input  and  loss  mechanisms  of  the  device.  Since  a  typical 
low  power  arcjet  nozzle  is  on  the  order  of  1  cm  in  length,  with  particle  speeds  between  S-10 
km/sec  the  charged  particle  residence  time  (~10^  sec)  is  much  smaller  than  the  recombination  time. 
As  a  result,  the  flow  is  in  chemical  nonequilibrium,  leading  to  frozen  flow  losses  and  a  decrease  in 
arcjet  performance.  Other  factors  that  limit  arcjet  performance  include  thermal  losses  to  the  anode 
due  to  ohmic  heating,  current  attachment,  convection  and  radiation  to  the  wall,^  wall  friction  losses 
in  the  nozzle  (due  to  large  viscous  and  thermal  boundary  layer  growth),  nonuniform  profile  losses, 
electrode  erosion  and  arcjet  instabilities.  The  nozzle  flow  is  very  viscous  (Re  <  1000)  and  highly 
nonequilibrium,  with  electron  temperatures  typically  much  greater  than  ion  and  neutral 
temperatures. 

To  achieve  high  speciHc  impulse  and  thrust  efficiency,  the  specific  energy  P/ih  (input 
electrical  power  divided  by  propellant  mass  flow  rate)  must  be  increased.  At  high  P/m  the  anode 
overheats,  so  that  improvement  of  arcjet  performance  for  all  power  levels  relies  on  understanding 
the  interaction  between  the  plasma  and  the  anode  wall.  ^ 

In  both  the  constrictor  and  nozzle  regions,  energy  transport  from  the  plasma  to  the  walls  is 
a  strong  function  of  the  electron  number  density  and  temperature  distributions.  Nonequilibrium 
processes  determine  the  spatial  variation  in  current  attachment  and  anode  heating.  Such  processes 
dominate  in  the  thick  boundary  layer,  so  that  measurements  of  local  rie  and  Te  inside  the  boundary 
layer  will  be  helpful  in  understanding  the  current  attachment  and  energy  deposition  between  the 
plasma  and  the  anode.  The  disuibutions  of  current  density,  electric  Held  and  voltage  are  also 
poorly  understood,  and  are  strongly  coupled  to  the  ne  and  Te  distributions  near  the  constrictor  and 
anode  walls. 

Due  to  the  thick  viscous  and  thermal  boundary  layers  present  in  the  nozzle,  the  plasma  core 
flow  and  boundary  layer  processes  are  coupled.  The  volumetric  ohmic  heating,  j2/o,  determines 
where  the  energy  is  deposited,  which  in  turn  determines  the  electrical  conductivity  o(ne,Te).  This 
electrical  conductivity  distribution  then  prescribes  the  voltage  and  current  distributions,  which 
affect  the  arc  attachment  and  the  thermal  losses  to  the  anode.  Also,  viscous  effects  in  the  boundary 
layer  are  coupled  with  the  plasma  core  flow. 

*Ph.D.  Candidate,  Nicholas  T.  Tiliakos,  leads  our  research  group  in  developing  internal  diagnostics  for  the  arcjet. 
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3 . 2  Research  Approach 


The  main  objective  of  this  research  is  to  understand  the  fundamental  physical  processes 
governing  energy  deposition  in  the  constrictor  and  nozzle  regions  of  a  low  power,  2  kW  hydrazine 
arcjet.  This  involves:  ( 1 )  using  Langmuir  probes  for  internal  arcjet  diagnostics,  (2)  determining  the 
ohmic  heating  distribution  in  the  arcjet  and  what  conditions  govern  its  distribution,  (3) 
understanding  the  physical  mechanisms  concerning  current  attachment  and  therefore  anode  heating, 
and  (4)  validating  the  UIUC  arcjet  model  (Section  1). 

To  accomplish  these  objectives  local  electron  number  density  and  electron  temperature 
measurements  using  flush-mounted  Langmuir  probes  will  be  made  at  several  axial  and  azimuthal 
locations  in  the  nozzle  and  constrictor  regions.  Current  density  measurements  will  also  be  made, 
while  using  the  Langmuir  probes  as  conducting  pins.  In  addition,  spectroscopic  measurements 
will  be  performed  in  the  nozzle  plasma  core  flow  (see  Section  3.5).  Flush-mounted  Langmuir 
probes  will  not  only  obtain  data  in  the  constrictor,  but  also  near  the  constrictor  exit,  1 .5  mm 
downstream.  Such  data  (rie,  Te.  j)  in  close  proximity  to  the  constrictor  is  very  important,  since  it 
has  been  found  that  the  arcjet  operation  is  stable  when  the  arc  attaches  in  the  supersonic  region  just 
downstream  of  the  constrictor  exit. 

3.3  Previous  Research 

The  first  experimental  investigations  of  the  nozzle  region  of  a  low  power  arcjet  were 
conducted  at  the  NASA-Lewis  Research  Center.  Zube  et.  al  2  were  the  first  to  perform  emission 
spectroscopy  of  the  plasma  flow  inside  the  nozzle  of  a  1  kW  hydrazine  arcjet.  Several  0.25  mm 
holes  were  drilled  into  the  diverging  section  of  the  nozzle,  providing  optical  access  to  the  internal 
flow.  Measurements  included  atomic  electron  excitation,  vibrational  and  rotational  temperatures 
for  the  expanding  plasma  using  the  relative  line  intensity  techniques.  However,  no  data  were 
obtained  inside  the  constrictor  region  and  the  closest  data  point  was  2.4  mm  downstream  of  the 
constrictor. 

To  date,  emission  spectroscopy  has  been  the  main  diagnostic  used  for  internal  probing  of 
the  arcjet.  Measurements  of  plasma  conditions  inside  the  constrictor  have  been  performed  on 
medium  power  arcjets  (5-20  kW).3.4  jjje  work  of  Zube  et.al  ^  the  constrictor  had  a  2  mm  radial 
hole  covered  by  a  quartz  window  for  optical  access.  Spectroscopic  measurements  included  radial 
profiles  of  H2  excitation  temperature,  electron  density  and  atomic  hydrogen  density  in  the 
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constrictor.  However,  due  to  problems  with  line  width  determination  near  the  walls  only  the  data 
from  the  iimer  60-70  %  of  the  constrictor  diameter  gave  reliable  results.^  A  problem  with  this  type 
of  configuration  is  heat  transfer  and  temperature  restrictions  on  the  window,  which  make  it 
unrealistic  to  relate  such  results  to  flight-type  thrusters.^ 

Curran  et.al  ^  have  studied  arc  energy  deposition  in  the  nozzle,  the  nature  of  arc  attachment, 
and  its  effects  on  performance  characteristics  of  the  device.  They  used  a  1-2  kW  arcjet  that 
incorporated  a  segmented  anode  nozzle  with  5  axial  conducting  segments  isolated  from  each  other 
by  boron  nitride  spacers.  The  current  to  each  segment  and  the  potential  difference  between  each 
segment  and  the  cathode  were  measured.  Tlie  current  was  found  to  attach  diffusely  to  the  nozzle 
wall,  with  more  than  50%  of  the  input  power  added  in  the  diverging  section  of  the  nozzle.  The 
current  distribution  was  found  to  be  dependent  on  the  mass  flow  rate.  We  will  perform  current 
density  measurements  similar  to  Curran,  using  conducting  pins  (~.01S"  diameter)  embedded  in  a 
monolithic  anode,  thus  minimizing  disruptions  to  the  ano^?  heat  flow. 

Very  little  experimental  work  has  focused  on  investigating  azimuthal  symmetry'  in  the 
current  density  distribution  inside  a  low  power  arcjet.  Harris  et.al  ^  have  studied  anode  arc  motion 
in  a  water-cooled  radially  segmented  30  kW  arcjet.  It  was  found  that  the  azimuthal  current 
distribution  in  an  arcjet  nozzle  is  asymmetric.  With  our  conducting  pin  configuration,  we  will  be 
able  to  measure  azimuthal  current  density,  thus  determining  whether  the  current  distribution  is 
symmetric  in  the  2  kW  arcjet. 

3.4  Langmuir  Probe  Theory 

Langmuir  probes  have  been  a  powerful  means  of  plasma  diagnostics  since  the  early 
1920's,  when  the  method  was  first  implemented  by  Irving  Langmuir  and  his  co-workers.^  This 
type  of  probe  consists  of  a  small  metallic  electrode  (wire)  inserted  into  a  plasma.  The  probe  is 
connected  to  a  variable  power  supply,  allowing  the  probe  to  be  biased  at  various  positive/negative 
voltages  with  respect  to  the  plasma.  The  current  collected  by  the  wire  is  then  measured  as  a 
function  of  the  biasing  probe  voltage,  generating  a  probe  characteristic  (Fig.  3.1),  and  providing 
information  on  the  plasma  conditions. 

The  electron  number  density  and  electron  temperature  can  be  obtained  from  a  probe 
characteristic,  assuming  collisionless,  thin  sheath,  Langmuir  probe  theory.  A  thin  sheath  means 
that  the  probe's  sheath  thickness  is  much  smaller  than  the  probe  radius.  The  probe's  sheath 
thickness  is  on  the  order  of  several  Debye  lengths,  Xp.  The  Debye  length  is  the  distance  over 
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which  the  perturbing  effects  of  a  charge  will  tend  to  penetrate  into  a  plasma.  For  ne=2xl020  m-3, 
and  Te=1.9  eV  (-22,000  °K),  X,d=.71  microns,  which,  compared  to  a  .015"  (381  microns)  probe 
diameter,  means  that  this  is  a  thin  probe  sheath. 

Region  'A',  (Fig.  3.1),  is  called  the  electron  saturation  region,  where  the  probe  voltage, 
Vp,  is  greater  than  the  plasma  space  potential.  Vs,  thereby  attracting  electrons  to  the  probe  and 
repelling  ions.  A  layer  of  charge  builds  up  on  the  probe  forming  an  electron  sheath.  In  the  sheath 
layer  the  plasma  is  no  longer  quasi-neutral,  the  number  density  of  electrons  and  ions  are  not  equal 
so  that  large  electric  fields  exist.  Outside  the  sheath  the  electric  field  is  small  so  that  the  plasma  is 
essentially  undisturbed  and  quasi-neutral.  Region  'A'  provides  the  magnitude  of  the  electron 
saturation  current.* 


Fig.  3. 1  Schematic  of  a  tj^ical  probe  characteristic,  showing  the  various  regions  of 
interest  that  provide  information  on  the  plasma  conditions.  [Reference  8] 


When  the  probe  voltage  equals  the  plasma  space  potential.  Vs,  there  is  no  net  voltage 
gradient  around  the  probe,  and  no  sheath  exists  that  repels  the  electrons.  As  a  result  the  charged 
particles  migrate  to  the  probe  because  of  their  random  thermal  velocities.  In  region  'B',  the 
Transition  Region,  the  probe  voltage  is  made  negative  relative  to  the  plasma  space  potential  so  that 
low  energy  electrons  are  repelled  and  ions  are  accelerated  to  the  probe.  If  the  electron  distribution 
were  Maxwellian,  then  the  shape  of  the  Transition  Region  on  a  log(Iprobc  -  Iion.sat.)  versus  Vp  plot 
would  be  linear,  with  the  slope  providing  the  electron  temperature,  as  given  in  Eq.  (3.1).®  Once 
the  electron  temperature  is  known,  the  electron  number  density  is  found  from  the  magnitude  of 
either  the  ion  or  electron  saturation  current,  as  given  in  Eq.  (3.2):® 


Tc  =  ^ 


I 


lon.sat 


(dl/dVo) 


.) 


(3.1) 


62 


ne  = 


I 


jgn.sati 


(0.61)eA^Tc/ini) 


1/2 


(3.2) 


where  Ap  is  the  probe  area,  lion.sat  is  the  ion  saturation  current,  mj  is  the  ion  mass  and  dl/dVo  is 
the  slope  of  the  probe  characteristic  of  the  log(Iprobc  -  lion.sat.)  versus  Vp  plot  (in  Region  'B'). 


When  the  probe  potential  equals  the  floating  potential,  Vf,  the  flux  of  electrons  equals  the 
flux  of  ions  to  the  probe  so  that  the  net  current  is  zero.  This  is  equivalent  to  inserting  a  wire, 
without  applying  a  biasing  voltage,  into  a  plasma. 


In  region  'C,  the  ion  saturation  region,  almost  all  electrons  are  repelled,  resulting  in  an  ion 
sheath.  The  ion  saturation  current  is  much  less  than  the  electron  saturation  current  because  the 
electron  mass  is  ~  1000  times  smaller  than  the  ion  mass  and  Te»Ti. 


plasma . 
potential 


potential  distribution 


x  =  0 


Fig.  3.2  Schematic  of  the  anode-plasma  interface  depicting  sheath  formation  and  the 
presheath  transition  region.  This  sheath  formation  is  for  a  grounded  anode, 
surrounded  by  an  electropositive  plasma. 


Figure  3.2  shows  the  sheath  formation  at  the  anode.  A  sheath  forms  at  the  anode-plasma 
interface  due  to  a  potential  drop  between  the  anode  and  the  plasma.  Near  an  electrode  there  are  not 
enough  collisions  to  randomize  the  electron  and  ion  flux,  and  a  sheath  forms  to  shield  the  bulk 
plasma  from  this  region  of  non-neutrality.  There  is  also  a  transition  region  between  the  sheath  and 
the  plasma  known  as  the  presheath.  In  order  to  satisfy  Poisson's  equation  in  the  sheath,  the  ions 
must  enter  the  presheath  with  the  Bohm  velocity,  Ui,Bohm-  The  potential  drop  at  the  anode/plasma 
interface  can  be  as  much  as  -15  Volts. 
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Langmuir  probes  are  easy  to  implement  since  they  are  relatively  simple  devices,  but 
interpreting  the  probe  characteristic  can  be  difficult,  especially  if  there  are  ionization  and  collisional 
reactions  occurring  inside  the  sheath.  If  the  probe  sheath  is  collisional  (i.e.  if  the  charged  particle 
mean  free  path  «  Xd)  then  the  Langmuir  probe  theory  requires  modification,  since  collisions 
reduce  the  current  to  the  probe  because  particles  diffuse  up  to  the  probe  instead  of  arriving  there  by 
free  flight.^ 

Figure  3.3  shows  a  schematic  of  plasma  flow  over  a  flush-mounted  Langmuir  probe  inside 
the  arcjet  nozzle.  There  are  three  regions  of  interest;  (1)  the  nozzle  viscous  boundary  layer,  (2)  the 
anode  sheath,  and  (3)  the  probe  sheath.  Initially,  when  the  probe  is  not  biased,  the  anode  and  the 
probe  sheaths  coincide.  As  the  probe  is  biased  positive  or  negative  a  separate  sheath  forms, 
surrounding  the  probe.  The  Langmuir  probe  measurements  of  ne  and  Te  are  thus  local 
measurements,  i.e.  inside  the  anode  boundary  layer. 


Fig.  3.3  This  figure  shows  the  plasma  flow  over  the  flush-mounted  Langmuir  probe 
inside  die  boundary  layer.  A  separate  anode  and  probe  sheath  form. 

A  variation  on  the  Langmuir  probe  is  the  double  Langmuir  probe,  which  consists  of  two 
probes  connected  to  a  biasing  voltage  supply  with  typical  voltage  sweep  rates  on  the  order  of  100 
jisec.lO  Analysis  of  the  characteristic,  similar  to  the  single  probe,  provides  measurements  of  ng 
and  Tg.^  The  benefit  of  a  double  Langmuir  probe  (or  a  "floating"  double  probe)  is  that  it  minimizes 
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electron  current  collection.  Consequently,  the  total  current  to  the  probes  can  never  be  greater  than 
the  ion  saturation  current,  minimizing  the  disturbance  on  the  plasma  discharge,  and  preventing 
oveibeating  of  the  probe. 

3.5  Experimental  Approach 


A  monolithic  anode  2  kW  hydrazine  arcjet  thruster  has  been  fabricated  at  the  University  of 
Illinois  (Fig.  3.4).  The  thruster  anode  is  machined  from  a  tungsten  alloy,  HD  18  (95%  W,  3.5% 
Ni,  1.5%  Fe)  which  is  easier  to  machine  than  thoriated  tungsten,  the  conventional  arcjet  anode 
material.  The  monolithic  anode  eliminates  the  front  seal,  simplifying  heat  transfer  analysis. 
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Fig.  3.4  The  new  UIUC  2  kW  hydrazine  arcjet  thruster. 

To  obtain  electron  number  density  and  electron  temperature  measurements  in  the  arcjet 
boundary  layer,  flush-mounted  single  and  double  Langmuir  probes  will  be  installed  in  the 
constrictor  and  nozzle  wall.  The  constrictor  diameter  has  been  increased  to  1.0  mm,  up  from  0.64 
mm  for  the  1  kW  thruster,  and  the  constrictor  length  is  now  1.5  mm,  up  from  .25  mm.  The  larger 
scale  thruster  can  incorporate  an  array  of  Langmuir  probes,  as  shown  in  Fig.  3.5.  The  probes  are 


65 


mounted  to  the  nozzle  through  vacuum  sealed  male  connectors.  The  2  kW  thruster  is  seated  on  a 
mounting  pod  inside  a  1.5  m^  vacuum  tank,  with  a  background  pressure  of  ~  100  mTorr  for  flow 
rates  of  SO  mg/sec. 


Fig.  3.S  Cross-sectional  schematic  of  a  typical  arcjet  configuration,  with  an  array  of 
Langmuir  probes  inside  the  constrictor  and  nozzle  walls. 

Once  thruster  operating  characteristics  have  b^n  determined,  the  anode  body  will  undergo 
additional  machining  to  accomodate  the  single  and  double  flush-mounted  Langmuir  probes.  Initial 
experimental  tests  will  implement  the  single  probes.  These  probes  are  made  of  tungsten  wire, 
surrounded  by  an  alumina  insulator,  and  inserted  in  a  stainless  steel  tube.  In  the  constrictor  there  is 
one  probe  with  diameter  .007  "  (.  18  mm),  while  m  the  nozzle  waU  there  can  be  at  most  six  probes, 
two  at  the  same  axial  location,  but  different  azimuthal  location  (separated  by  90°),  with  diameter 
.015  "  (0.38  mm).  The  probes  will  be  biased  with  a  function  generator  sweeping  between  ±12 
volts.  A  power  op-amp  will  be  be  used  to  isolate  the  probe  from  the  function  generator  in  order  to 
avoid  drawing  excessive  current  from  the  function  generator. 

The  next  phase  of  the  experiment  will  involve  using  double  Langmuir  probes  and  perhaps 
implementing  more  probes  at  different  locations  inside  the  nozzle.  To  obtain  the  distribution  of 
ohmic  heating  in  the  nozzle  region  the  current  density  will  be  measured  along  the  nozzle  wall.  This 
is  done  by  measuring  the  current  to  each  of  the  Langmuir  probes,  which  for  this  case  are  now  just 
conducting  pins  connected  directly  to  the  anode.  Knowing  the  cross-sectional  area  of  each  pin  the 
current  density  distribution  along  the  nozzle  is  obtained.  The  current  density  measurements, 
coupled  with  calculations  of  the  electrical  conductivity  distribution  (based  on  ne,  Te  measurements) 
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provide  an  estimate  of  the  ohmic  heating,  j2/o(nc,  Te).  at  various  axial  locations  inside  the 
boundary  layer.  Calculations  of  the  electrical  conductivity  will  also  lead  to  the  voltage  distribution 
along  the  nozzle  wall. 

Spectroscopic  measurements  may  eventually  be  required  to  obtain  ng  and  Te  profiles  in  the 
plasma  core  flow  in  the  nozzle.  These  measurements  would  be  made  through  a  series  of  open 
holes  in  the  nozzle  wall.  Emission  spectroscopy  will  require  use  of  an  Optical  Multichannel 
Analyzer  (OMA).  Other  diagnostics  would  include  static  pressure  measurements  along  the 
constrictor  and  nozzle  walls. 

The  final  phase  of  the  experiment  will  require  modifying  the  UIUC  arcjet  model  to  include 
applied  magnetic  field  effects  and  observing  how  the  field  affects  n^,  Tc  and  j  distributions.  Based 
on  the  numerical  results  a  magnetic  bobbin  can  be  designed  and  used  as  a  diagnostic  tool. 

Lai  t  year  an  investigation  of  a  pulsed  magnetic  field  on  the  operation  of  a  NASA  1  kW 
arcjet  was  :onducted  in  the  1.5  m^  tank  at  the  UIUC  Electric  Propulsion  Laboratory.  ^  We  found 
that  during  the  duration  of  the  pulse  the  arcjet  voltage  increased  by  as  much  as  6  volts.  Similarly, 
for  a  given  magnetic  field  strength,  as  the  propellant  flow  rate  was  decreased  the  change  in  arcjet 
operating  volmge  due  to  the  magnetic  field  increased.  This  preliminary  work  suggested  that  a 
magnetic  field  could  be  used  to  control  the  current  attachment  location,  i.e.  the  electron  drift 
motion,  and  therefom  the  he  ating  inside  the  arcjet. 

Some  possible  effects  of  an  applied  magnetic  field  include  decreasing  thermal  losses  to  the 
anode  wall,  since  the  electron  thermal  conductivity  scales  as  fV(l  +  where  is  the  ele''tron 
Hall  parameter  (proportional  to  the  magnetic  field  strength).  Also,  propellant  swirl  would  be 
increased  due  to  the  jrBj  and  jzBr  Lorentz  forces;  this  would  enhance  mixing  and  arc  stability. 
Similarly,  particle  residence  time  would  increase,  allowing  more  time  for  recombination  to  occur.  •  ^ 

A  steady  state  applied  field,  using  a  pulsed  magnet  in  series  with  the  arcjet,  can  be  used  to 
modify  the  current  distribution  in  a  predictable  way,  and  can  therefore  assist  in  further  validating 
the  UIUC  arcjet  model  (Section  1).  Magnetic  field  strengths  of  several  100  Gauss  can  be  achieved 
on  the  nozzle  centerline.  The  UIUC  code  will  be  modified  to  include  the  effects  of  an  externally 
applied  magnetic  field.  This  will  involve  incorporating  analytic  expressions  for  the  externally 
applied  magnetic  field,  which  has  components  Br  and  Bj. 

This  investigation  will  involve  use  of  the  UIUC  arcjet  mo..  comparisons  with  the 
experimental  results.  The  numerical  model  will  also  assist  in  directing  the  experimental  work. 
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Initially,  the  model  will  be  used  to  predict  electron  number  density,  electron  temperature  and 
current  density  distributions  in  the  nozzle  region  for  various  propellant  flow  rates  and  propellant 
mixtures.  The  mole  fraction  of  N2,  x,  is  also  an  input  parameter,  since  the  propellant  can  be 
described  as  XN2  +  H2.  If  x  =0  then  the  propellant  is  H2,  while  if  x  =0.5  the  propellant  is 
simulated  hydrazine.  This  leads  to  a  flexible  model  that  can  predict  arcjet  performance  for 
ammonia,  hydrogen  and  simulated  hydrazine  propellants.  This  variation  in  propellant  type  can  be 
implemented  in  the  laboratory  as  well.  The  model  predictions  for  rie,  Te,  j  and  a  distributions  will 
then  be  compared  with  the  experimental  data  in  the  constrictor  and  nozzle  regions. 

3.5.1  Experimental  Procedure 

Figure  3.6  shows  a  schematic  of  the  probe  experiment.  The  probe  supply  consists  of  a 
function  generator  that  sweeps  the  Langmuir  probe  between  ±15  V.  A  power  op-amp  will  be  used 
as  an  isolation  stage  between  the  probe  and  tl^  function  generator.  The  tests  wUl  be  run  during  the 
20  minute  arcjet  warm-up  period  to  minimize  probe  thermal  loading. 


Fig.  3.6  Schematic  of  Langmuir  probe/conducting  pin  experiment. 

The  probe  current  is  obtained  by  measuring  the  voltage  drop  across  a  known  resistance, 
while  the  probe  voltage  is  measured  by  a  voltage  probe  connected  to  a  10  mHz  SOLTEC 
oscilloscope.  The  SOLTEC  oscilloscope  will  be  set  to  acquire  data  for  a  period  at  least  as  long  as  it 
takes  the  function  generator  to  sweep  ±15V.  The  data  it  will  obtain  is  probe  voltage  vs  time  and 
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AVr  vs  tir  where  AVr  is  the  voltage  drop  across  a  known  resistance,  -  probe  current.  The  two 
channels  tor  probe  voltage  and  AVr  vs  time  are  cross  plotted  to  obtain  the  probe  characteristic. 

3.5.2  Some  Important  Probe  Considerations 

Some  of  the  consequences  of  physical  probing  of  the  hostile  arcjet  environment  include 
probe  contamination,  probe  surface  reactions  and  probe  heating. 

PrQt)g  CQntamination/Syrf&gs  CQnsideraUons 

( 1 )  Thermionic  emission  of  electrons  from  the  probe's  surface  can  lead  to  a  perturbation 
to  the  probe’s  electric  field.  As  a  result  these  emission  currents  could  greatly  perturb  the  potential 
distribution  in  the  neighborhood  of  the  probe,  thereby  invalidating  any  probe  theories  used  that  do 
not  take  emission  currents  into  account.  In  order  for  thermionic  emission  effects  to  be  negligible 
the  probe  temperature  should  never  exceed  1200  °C  (1473  '’K).l2  This  may  be  particularly 
important  for  probes  near  the  constrictor  region  where  the  temperatures  are  expected  to  be  highest. 

(2)  If  the  probe  is  at  a  different  temperature  than  the  plasma  this  could  lead  to  a 
perturbation  to  the  concentration  and  temperature  of  the  carriers  and  neutrals.  As  a  result  of  these 
temperature  gradients,  heat  and  mass  flow  occurs.  This  should  not  be  a  problem  for  our  design 
since  the  probe  is  insulated  in  alumina,  which  conducts  heat  well,  and  also  because  the  probe  is 
small  relative  to  the  nozzle  wall  in  which  it  is  embedded;  as  a  result,  the  probe  will  approximately 
equal  the  anode  wall  temperature.  The  nozzle  wall,  in  turn,  is  in  good  thermal  contact  with  the 
plasma. 


Since  the  probe  tip  receives  the  greatest  energy  flux  it  is  usually  made  out  of  a  refractory 
metal  such  as  platinum,  molybdenum  or  tungsten.^  In  our  investigation  the  probe  will  be  made  of 
tungsten. 

Problems  with  the  probe's  surface  can  lead  to  plasma  perturbations  resulting  in  errors  in  the 
measurement  of  current  density  or  probe  potential.  Some  of  these  problems  include: 

(1)  The  probe  surface  may  be  involved  in  direct  chemical  reaction  or  catalysis  in 
chemically  active  plasmas,  resulting  in  perturbations  in  both  plasma  composition  and  temperature. 
This  is  not  a  problem  for  us  since  the  plasma  is  not  chemically  reacting  with  the  probe  surface. 
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(2)  Sputtering  of  the  probe's  surface  material  can  lead  to  uncertainties  in  the  probe's 
effective  collecting  area,  carrier  reflection  and  emission  from  the  surface,  leading  to  errors  in  the 
current  density  measurement.  If  the  sputtered  material  is  in  electrical  contact  with  the  probe  there 
will  be  an  under-estimation  of  the  effective  collecting  area  of  the  probe.  Similarly,  if  the  insulating 
material  sputters  onto  the  probe's  surface  there  could  be  an  over-estimation  of  the  effective 
collecting  area.  One  way  to  avoid  problems  with  changing  the  effective  area  of  the  probe  is  to 
center  the  probe  tip  in  the  insulator  so  that  the  point  of  contact  with  the  insulator  is  well  away  from 
the  plasma.^ 


Probe  Heating 


In  considering  sizing  of  the  probe  and  the  experimental  procedure  of  obtaining  the  probe 
characteristic,  we  must  account  for  the  maximum  heat  load  that  the  probe  can  tolerate.  This  is 
particularly  important  if  the  electron  saturation  portion  of  the  characteristic  is  obtained. 

To  avoid  over-heating  the  probe,  especially  when  collecting  electrons,  the  probe  voltage 
will  be  pulsed.  The  pulse  begins  at  t  =  0,  resulting  in  an  increase  in  probe  temperature,  AT^, 
given  by  subscript  'w',  which  refers  to  tungsten  probe  properties.  A  simple  calculation  of  the 
temperature  rise  profile  of  the  probe  is  given  by  the  following  equation:*^ 


ATw  = 


2qf^ 

(iCPvyCwKw) 


(3.3) 


where  q  is  the  heat  flux  to  the  probe  assumed  to  be  constant  and  equal  to: 


q  = 


_  ^ e,satfe,sat  _  ^ p^e.sat 
Aprobc  Aprobe 


[W/m2] 


(3.4) 


It  is  assumed  that  the  electron  saturation  voltage  is  approximately  equal  to  the  plasma 
potential,  which  is  very  small  (~  0(1)  mV).  The  calculation  is  made  for  conditions  just 
downstream  of  the  constrictor  and  for  the  electron-saturation  region  of  the  probe  characteristic, 
which  are  worst  case  conditions.  The  results  are  shown  in  a  plot  of  ATprobe  versus  time  for  three 
different  plasma  potentials.  Fig.  3.7.  The  results  suggest  that  we  should  pulse  the  probe  on  and 
off  when  near  the  electron  saturation  region,  because  the  probe  may  heat  up  to  as  much  as  1250 
°K  -I-  ATprobe  or  about  1500  °K.  The  probe  must  never  exceed  about  1473  "K  (1200  °C)  to  ensure 
that  thermionic  emission  from  the  probe  is  negligible. 
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Temperature  Rise  of  Probe  [‘K  ] 


Temperature  Rise  Profile  of  Single  Probe 

[Near  Constrictor;  n.«2.1049xl0“  m-*,T^-22,000  “K  ] 
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Fig.  3.7  Temperature  rise  profile  of  a  single  Langmuir  probe  for  various  assumed 
plasma  potential  values.  To  prevent  excessive  probe  heating  the  probe 
biasing  voltage  should  be  swept  at  a  frequency  corresponding  to  ATw=:20  °. 
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Appendix  A  Nitrogen  Ion  Density  in  Non-Equilibrium  Hydrazine  Arcjets 

It  has  been  observed  experimentally  by  Keefer  et  al.^  and  Hargus  et  al.^  that  unusually  low 
nitrogen  densities  are  found  downstream  of  the  constrictor  in  ammonia  and  hydrazine  arcjets.  It  is 
suggested  that  this  is  a  result  of  charge  exchange  collisions  in  the  expanding  nitrogen-hydrogen 
propellant. 

Number  density  ratios  for  nitrogen  and  hydrogen  atomic  and  ionic  species  are  shown  in 
Fig.  A.l.  Calculations  are  for  a  weakly  ionized  hydrazine  plasma  in  thermal  (Te  =  Tg)  and 
chemical  equilibrium  at  p  =  latm.  For  gas  temperatures  below  Tg  =  5000,  the  ion  ratios  are  less 
thanN+/H+=  0.05.  This  is  due  to  the  preferential  dissociation  of  H2  (Ed  =  4.48  eV)  over  N2  (Ed 
=  9.61  eV)  indicated  in  Fig.  A.2,  where  molecular  and  atomic  number  densities  are  given.  The 
relatively  low  nitrogen  ion  densities  observed  near  the  nozzle  exit  suggest  that  the  effective  frozen 
gas  temperature  is  less  than  Tg  =  6000  K. 


Simulated  hydrazine  at  p  s  1  atm 


Temperature,  K 


Figure  A.1:  Ionic  and  atomic  number  density  ratios  for  simulated  N2H4  at  p  =:  1  atm. 
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Figure  A.2:  Molecular  and  atomic  number  densities  for  simulated  N2H4  at  p  =  1  atm. 

The  arcjet  plasma  is  modeled  as  a  highly  dissociated,  partially  ionized,  low-pressure  N-H 
plasma,  with  Te  >  Tg  and  a  Tg  of  4000  <  Tg  <  10,000  K  in  steady  state.  Under  these  conditions 
the  hydrogen  is  more  highly  dissociated  than  the  nitrogen  and  fnxn  equilibrium  calculations  nN  is 
negligible  for  Tg  <  6000  K.  Furthermore,  ionization  by  electron  collision  is  an  infrequent  process 

Under  these  conditions,  nitrogen  ions  r^idly  lose  their  charge  by  charge-exchange 
collisions  with  hydrogen,  the  cross-section  for  which  is  on  the  same  order  of  magnitude  as  N-H 
momentum  transfer  collisions.  Thus; 


H  +  N+  — >  H+  +  N  and  N  +  H+  — >  N+  +  H 


Since  the  charge  exchange  collision  rate  is  » the  electron  ionization  rate,  these  two  reactions 
balance  in  steady  state. 


74 


Writing  the  volumetric  charge  exchange  collision  rate  in  terms  of  a  charge  exchange  cross 
sections  Qh-n+  and  Qn-h+>  with  a  relative  thermal  speed  vh,n=vh: 

[collisions/m^-s]  =  nnnN+Qn-N+VH  =  ONnH+QN-H+VH 

Assuming  that  Qh-n+  =  Qn-h+  and  combining: 

Pn*  _  Hn 
Ph"  PH 

For  hydrazine  arcjet  plasmas,  n^  «  n^,  implying  a  low  nitrogen  ion  density.  For  plasmas 
with  low  gas  temperature  Tg  «  Tg,  and  Te  <  10,000  K,  the  nitrogen  ion  density  is  effectively 
determined  by  the  gas  temperature  Tg  and  not  the  electron  temperature.  Thus  as  Tg  falls  rapidly  in 
the  nozzle  by  expansion  and  heat  loss,  the  nitrogen  ion  density  nN+  also  falls,  reaching  low  values 
at  the  exit  plane. 
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